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SUMMARY 


The NASA Langley Research Center has been conducting research to establish founda- 
tions for adequate representation and treatment of the airframe structure in design 
analyses of helicopter vibrations. This paper presents a body of formulations to 
couple a finite-element analysis model of the airframe to a rotor analysis model and 
calculate airframe vibration levels. The rotor is represented by a general set of 
linear differential equations with periodic coefficients, and the connections between 
the rotor and airframe are specified through general linear equations of constraint. 
Background is provided relating to development of rotor and airframe models, and as an 
aid to structural engineers, the origins of linearized rotor equations are reviewed. 
Coupling equations are derived and then applied to combine the rotor and airframe 
equations into one set of linear differential equations governing vibrations of the 
rotor-airframe system. These equations are solved by the harmonic balance method to 
yield the system steady-state vibrations. A key feature of the solution process is to 
represent the airframe in terms of forced responses calculated at harmonics of the 
rotor rotational frequency. A method based on matrix partitioning is presented for 
quick recalculations of vibrations in design studies when only relatively few airframe 
members are varied. All relations are presented in forms suitable for direct computer 
implementation. During a discussion of the method, it is argued that a properly lin- 
earized rotor model should be as good for coupled rotor-airframe vibration analysis 
as the underlying nonlinear model, and representing the airframe by harmonic forced 
responses is recommended for airframe structural design work. 


INTRODUCTION 

Helicopters are prone to vibrations arising from the cyclic nature of the rotor 
actions. The vibrations normally pervade both the rotor and the airframe and can 
seriously degrade both service life and ride qualities. Vibrations also frequently 
limit the maximum speed in forward flight. In U.S. civil and military helicopters, 
it is usual to incorporate devices solely to attenuate vibrations. These devices have 
added significant weight and complexity. 

U.S. helicopter companies have rarely relied on analysis during design in their 
efforts to limit vibration. With only a few exceptions, helicopters have been 
designed to performance requirements with past experience with vibrations taken into 
account, and new problems with vibration have been solved during flight tests and 
operation. There is now a recognized need to account for vibrations during the 
analytical phases of design. The advent of modern methods of computer analysis has 
provided the opportunity to achieve such a capability. 

Four technical factors influencing vibrations of a helicopter should be recog- 
nized: (1) vibratory loads induced by the rotor actions, (2) response of the rotor, 

(3) coupling of the rotor and airframe, and (4) response of the airframe. The great- 
est excitation of vibrations is caused by cyclic loads on the rotor blades due to 
their interaction with the airstream. The dynamic characteristics of the rotor and 
airframe and the coupling of these two systems determine the manner in which the 
helicopter responds to this excitation. Thus, approaches to establish reliable 
methods for vibration analysis should address all four factors. 


Among the analysis methods now employed by industry applicable to helicopter 
vibrations, two categories can be distinguished: (1) methods for analysis of airframe 

behavior and (2) methods for analysis of rotor behavior. For the nonrotating airframe 
components, the NASTRAN® computer code (ref. 1) , which embodies the finite-element 
method for structural analysis, has become a standard tool used throughout the heli- 
copter industry for design calculations of internal structural loads and checks on 
airframe vibrations . For rotating components , there has been extensive work on for- 
mulating equations of motion of rotors and devising computer solutions of the equa- 
tions. (See, for example, refs. 2 to 6.) These categories of methods have evolved 
during the past 25 years since high-speed digital computers were introduced. The 
methods for rotors and airframes have generally been developed separately. To realize 
the ideal of a trusted vibration analysis, continued work is required to improve air- 
frame analysis methods and rotor analysis methods separately. In addition, there is 
a requirement for general methods for combining mathematical models of the rotor and 
the airframe to calculate the vibrations of a helicopter as a single system. 

The literature on practical methods for calculating vibrations of a helicopter 
as a single system has been sparse until quite recently when a number of papers have 
appeared. The earlier papers come from Gerstenberger and Wood (ref. 7) and Novak 
(ref. 8) . In their analyses, a key feature is representation of the airframe by 
impedances at the rotor attachment points. In a later paper, Staley and Sciarra 
(ref. 9) also applied this approach. Reference 10 describes a vibration analysis 
of a helicopter design notable for the exposure of the treatment of the airframe. 

Here the airframe is also represented by impedances. The work of these earlier 
investigators includes correlations of analytical results with measured flight vibra- 
tion data. References 7 to 10 provide a technical basis on which to formulate a 
general method of vibration analysis suitable for airframe structural design work. 
Three other recent papers (refs. 11 to 13) have specifically addressed practical 
methods for calculating helicopter vibrations. In contrast to references 7 to 10, 
they use a modal representation of the airframe, which is a significant alternative 
line of development. 

It is recognized that this discussion excludes a number of existing computer 
simulations of the helicopter in flight (see, for example, refs. 3, 4, 6, and 14). 

Such simulations, of course, incorporate representations of both the rotor and the 
airframe and the connections between the two and thus theoretically could be applied 
to calculate vibrations. However, there is little note in the literature of their 
use to calculate airframe vibrations. These simulations have been applied mainly to 
development of flight controls, to checks on the stability of rotors, and to calcula- 
tion of blade vibratory loads. As a rule, the current simulations incorporate only 
cursory, if any, treatment of the airframe elasticity, and they would be cumbersome 
to use directly for airframe structural design work. 

References 15 to 19 contain relevant subsidiary analysis procedures and solutions 
of vibrations of simplified rotor-airframe systems. Reference 15 is particularly 
interesting since it is evidently the first paper to pose the problem of generalizing 
the harmonic balance method and is evidently the only attempt to solve this problem 
prior to the present paper. The harmonic balance method has been widely used to solve 
for steady-state vibrations of particular rotating systems, but a perfectly general 
algorithm is needed to systematize analysis of rotor-airframe vibrations. Also, as a 
matter of interest, the reader may compare the present paper with reference 16 with 
regard to the assumptions underlying the equations which express the coupling of the 
rotor to the airframe . 
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The present paper is an outcome of recent efforts at the NASA Langley Research 
Center to establish foundations for adequate representation and treatment of the air- 
frame structure in design analysis of helicopter vibrations. The subject addressed 
is the problem of coupling a general finite-element model of the airframe system to 
rotor mathematical models. Presented are a body of formulations and discussions 
intended as a basis of both engineering and computational theory for efficiently 
implementing such couplings and calculating the airframe vibration levels. Emphasis 
is placed on analysis of vibrations during industrial design of airframe components. 
Rotor analysis is discussed, but with the limited purpose of ensuring correct inclu- 
sion of rotor effects in the airframe analysis. 

The presentation consists of a text and a number of appendices. The text 
describes the salient features of the proposed method. The appendices provide back- 
ground on equations which appear in the text and discuss subsidiary topics. The text 
covers the following topics: 

1. The general form of the equations for the three components of a linear 

coupled rotor-airframe vibration model: linearized rotor model, finite- 

element airframe model, and linearized interface model 

2. The steps for combining the component equations to form the coupled equations 

3. The steps for solution of the coupled equations by the harmonic balance 

method, including reduction of the number of airframe degrees of freedom in 
the harmonic balance equations by use of airframe harmonic forced responses 

4. A block diagram indicating the basic sequence of tasks for calculating air- 

frame vibrations 

5. The plausibility and utility of a linear rotor-airframe model 

6. Computational advantages of the method 
The appendices discuss the following: 

1. The three component models 

2. An alternative representation of the rotor by impedances 

3. Subsidiary computations 

4. A method for quick recalculation of airframe vibrations in design studies 

when only relatively few airframe members are varied 

5. Additional block diagrams indicating the sequence of tasks for calculating 

airframe vibrations including steps for reanalysis and for representation 
of the rotor by impedances 

PREPARATORY REMARKS 
The Helicopter Vibration Problem 

The system treated is a flexible helicopter in steady flight as sketched in 
figure 1. The vehicle is viewed as stationary and the air as moving. For convenience 
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Figure 1.- Steady-state vibrations of a flexible helicopter 

in steady flight. 


of discussion, it is assumed that there is only one rotor and that the required anti- 
torque is provided by constant external forces. A turning rotor interacting with the 
airstream gives rise to vibrations which normally pervade both the rotor and the air- 
frame. In steady flight, the vibration amplitudes take on a periodic character, as 
indicated in the figure. 


Steady Flight 

Examples of steady flight conditions are hover, straight and level flight, and 
circling with constant speed, bank angle, and altitude. Nonsteady conditions occur- 
ring during maneuvers such as vertical ascents, pull-ups, and push-overs can be 
realistically treated under steady flight assumptions if, as usual, the time to 
execute the maneuver encompasses many rotor revolutions. 


Rotor and Airframe Systems 

Figure 2 illustrates the helicopter components which must be considered in order 
to represent the interface between the rotating and nonrotating parts. The portion 
of the helicopter containing all those parts which in the analysis are considered to 
be rotating is called the rotor system. The remaining portion of the helicopter is 
called the airframe system. 
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Figure 2.- Components of a helicopter which typically must be considered 
in defining the interface between the rotor system and the airframe 
system. 


Figure 3 illustrates two ways in which rotor systems could be defined for analy- 
sis. In figure 3(a), the rotor system is limited to blades and hub. In this case, 
some parts of the airframe system must be considered to be nonrotating even though 
certain of these parts such as the rotor shaft, pitch links, and rotating swashplate 
actually do rotate. Thus, this point of view implicitly incorporates a fictitious 
bearing connecting the hypothetically stationary rotor shaft and the rotating hub. 

The torque needed to maintain the rotational speed of the rotor is implicitly assumed 
to be externally applied at the hub bearing. Effects in the rotor system resulting 
from rotation of components not included in the rotor system can be accounted for by 
specifying an impedance at the hub. A characteristic of the point of view of fig- 
ure 3 (a) is that only one load path leading from the rotor system to the airframe 
system, namely, that provided by the shaft, is taken into account. This point of view 
of the rotor system has been the one predominantly adopted in analysis of helicopter 
vibrations . 

In figure 3(b), the rotor system includes the blades, hub, rotor shaft, pitch 
links, rotating swashplate, nonrotating swashplate, control actuators, and trans- 
mission. The engine and engine shaft, which are included in the airframe system, are 
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Blades 


(a) System consisting of blades and hub. 



(b) System consisting of blades and hub, rotor shaft, pitch 
links, rotating swashplate, nonrotating swashplate, 
control actuators, and transmission. 

Figure 3.- Two illustrative definitions of the rotor system. 


considered to be nonrotating. Thus, the point of view of figure 3(b) implicitly 
incorporates a fictitious bearing in the engine shaft. The torque needed to maintain 
the rotational speed of the rotor is implicitly assumed to be externally applied at 
this fictitious bearing. The impedance effects of rotating components in the engine 
which are not included in the rotor system can be accounted for as an impedance at 
the fictitious bearing. The point of view of figure 3(b) allows additional load paths 
between the rotor and the airframe. 

In this paper, the restriction is made that the rotor system may include only 
those rotating parts having rotational speeds which are integer multiples of the rotor 
speed. This would generally rule out direct inclusion of the engine and engine shaft 
(fig. 2) as rotating parts, since they usually operate at a noninteger multiple of 
the rotor speed. Effects of engine rotation can be indirectly included as an imped- 
ance, as mentioned. 
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EQUATIONS OF MOTION OF COUPLED ROTOR-AIRFRAME SYSTEM 


Three elements are needed for forming the equations of motion of a coupled rotor- 
airframe system: (1) the rotor equations of motion, (2) the airframe equations of 

motion, and (3) the coupling equations. General linearized equations for these three 
elements are derived in appendices A, B, and C. In this section, these equations are 
briefly described, and then linear equations of motion are derived for a general 
coupled rotor-airframe system. 


Rotor Equations of Motion 

The linearized rotor equations derived in appendix A are presented in virtual 
work form as 



<Sw r dt 



6{dvr} t |{qro} + [mr]{dvr} 


+ [cr]{dVR> + [kr]{dVR}| dt = 0 


( 1 ) 


Equation (1) is a variational equation equivalent to a set of linear ordinary differ- 
ential equations of second order. The independent variable is the time t. In equa- 
tion (1) the matrices [mr] , [cr] , [kr], and {QRO} are known matrices which character- 
ize the rotor system. The elements of the vector {dvr} are unknown variables which 
determine the vibratory motions of the rotor. These variables are often called 
generalized coordinates. However, for visualization they may be thought of as repre- 
senting small displacements from a trim solution. The known matrices are all periodic 
functions of time with period 2TT/fi where £2 is the rotor rotational frequency. The 
matrices [mr], [cr], and [kr] represent the inertial, damping, and stiffness proper- 
ties of the rotor including contributions from both mechanical and aerodynamic models. 
The vector {QRO} comes from the rotor trim solution and is the vector of external 
forces required to maintain the conditions assumed in calculating the trim solution. 
Fourier series forms for the known matrices in equation (1) are defined as follows: 


[MR] = [MROC] + [MR1S] sin fit + [MRlc] cos fit + [MR2S] sin 2fit + ... (2a) 
[cr] = [CROC] + [CRls] sin fit + [CRlc] cos fit + [CR2s] sin 2fit + ... (2b) 
[kr] = [KROC] + [KRIS] sin fit + [KRlc] cos fit + [KR2s] sin 2fit + . . . (2c) 
{QRO} = {QROOC} + {QROIS} sin fit + {QROlc} cos fit + {QR02S} sin 2fit + ... (2d) 


The coefficient matrices appearing in these series, such as [MRls] and {QROls}, are 
constant . 

As discussed in appendix A, equation (1) is partitioned for use in coupling the 
rotor to the airframe : 
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6w dt 
r 



[MRll][MR12] 

[MR2l][MR22] 



[CRll][CR12] 

[CR21][CR22] 


|{dvb}| 
I (dvh}| 


+ 


[KRll][KR12] 

[KR2l][KR22] 



dt = 0 


(3) 


The submatrix {dVh} contains the rotor variables which appear explicitly in relations 
expressing connections of the rotor to the airframe. Of course, each of the sub- 
matrices appearing in equation (3) has a Fourier series representation analogous to 
equations (2). For example, 

[MRll] = [MROCll] + [MRlSll] sin St + [MRICll] cos St + ... (4a) 

{QBO} = {QBOOC} + {QBOls} sin St + {qBOIC} cos St + ... (4b) 


{QHO} = {QHOOC} + {OHOls} sin St + {QHOlc} cos St + ... 


(4c) 


In this paper, it is assumed that the linearized rotor equations of motion are 
specified at the outset of computation. The specification is made by 

1. Providing the coefficient matrices in equations (2) 

2. Providing the rotor rotational frequency S 

3. Identifying the variables appearing in {DVh} so that the partitions indicated 

in equation (3) can be made. 


Airframe Equations of Motion 

As discussed in appendix B, the airframe equations of motion are based on a 
typical airframe finite-element model. Formulated in terms of virtual work, the 
equations take the form, 



Sw„ dt = 


f 5{5}T { [ MA 




[ca]|z} + [ka]{z} + 



(5) 


The elements of the vector { z} are unknown variables which determine the 
of the airframe finite-element model. The matrices [ma], [ca] , and [ka] 


displacements 
represent the 
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mass, damping, and stiffness matrices of the model. These are known constant 
matrices. The vector {l} in equation (5) represents the loads directly applied to 
the airframe from sources which are external to both the airframe system and the rotor 
mechanical system. It is emphasized that {£} does not include loads imposed upon the 
airframe from the rotor system through mechanical connections . These rotor loads are 
treated as internal loads. Types of practical loads which can contribute to {l} 
include 


1. Constant aerodynamic and gravitational loads 

2. Oscillatory aerodynamic loads generated by the interaction of the rotor 

system with the airstream and impinging directly on the airframe 

3. Oscillatory test loads 

The following Fourier series form encompassing these types of loads is assumed 
for {l}: 


{l} = {LOC} + {LIS} sin fit + {Lie} cos fit + ... 


( 6 ) 


The coefficient vectors in this series are constant. These vectors, with the excep- 
tion of {LOC}, represent forces oscillating about zero mean and are assumed to be pre- 
scribed at the outset of computation. It turns out that it is not necessary to know 
the distribution of the mean loads acting on the airframe to calculate the steady- 
state vibration levels. However, it may be convenient for some purposes to designate 
static forc es acting on the airframe. In this paper, such known forces are repre- 
sent ed by {LOC}. Note that there is no requirement for the static forces represented 
by {LOC} to balance the mean rotor forces represented in appendix A by {qhtr}. 
Equilibrium of the mean rotor forces is explained subsequently. 

As noted in appendix B, a partitioned form of equation (5) is introduced to 
facilitate coupling the airframe to the rotor: 



dt 
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[MAll][MA12] 

[MA2l][MA22] 



[CAll][CA12] 

[CA2l][CA22] 



+ 


[KAll][KA12] 

[KA2l][KA22] 



dt 


0 


(7) 


In equation (7) , the vector {zh} contains the airframe variables which appear 
explicitly in relations expressing connections of the rotor system to the airframe 
system. The following Fourier series forms are specified for the partitions of {£} 
appearing in equation (7) : 
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{lh} = {LHOC} + {LH1S } sin fit + {LHIC} cos fit + ... 


(8a) 


{ls} = {LSOC} + {LSls} sin fit + {LSlc} cos fit + ... 


(8b) 


Coupling Equations 

The general coupling equations developed in appendix C are 


{dvh} = [th]{zh> 

(9a) 

[tc]{zh} = {o} 

(9b) 


where [th] and [tc] are known constant matrices. As noted in appendix C, the rows in 
both [th] and [tc] are linearly independent. The first equation explicitly expresses 
the rotor interface variables {dvh} in terms of the airframe interface variables, 
{zh}. The second equation expresses any constraints on the airframe interface vari- 
ables implied by the assumptions made in modeling the rotor system. If no such con- 
straints are implied, then equation (9b) is absent. 


Application of Constraints to Airframe Interface Variables 

To derive the coupled equations, the constraints represented by equation (9b) 
must first be applied to the airframe interface variables {zh}. This modifies both 
the first coupling equation (9a) and the airframe equation (7) and eliminates 
equation (9b) . 

Modification of first coupling equation .- When constraints are applied, some of 

the variables in {zh} are no longer independent variables. To distinguish dependent 
and independent variables, equation (9b) is written in partitioned form as 



{ 0 } 


( 10 ) 


where [tCd] is assumed to be a nonsingular square matrix. The vector {zhd}, which is 
associated with the submatrix [tCd] , represents the dependent variables. The vector 
{ZHl}, which is associated with the submatrix [TCl] , represents the independent vari- 
ables. Since the rows of [tc] are linearly independent, at least one nonsingular 
submatrix always exists. However, in general, elements of the vector {zh} must be 
renumbered and the columns of the matrix [tc] must be reordered to position such a 
submatrix to the left as has been done in equation (10) . 
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Equation (10) yields the following equation for elimination of the dependent 
variables : 



-[tcd] 1 [tci] 

[l] 


{ZHl} 


(ID 


The right side of equation (9a) is renumbered and reordered as was done in obtaining 
equation (10) from equation (9b) : 


{dvh} = 





( 12 ) 


Substitution of equation (11) into equation (12) yields the modified coupling 
equation 


{dvh} = [th]{zhi} 


(13) 


where 



[thd][tcd] 1 [tci]J 


(14) 


Modification of airframe equations .- As. discussed in appendix B, the application 
of linear constraints to finite-element models is straightforward and in the NASTRAN 
finite-element code has been rendered systematic. Nevertheless, the formal steps for 
incorporating constraints in a finite-element model are stated here. 

Equation (11) can be rewritten by reordering rows as follows: 


{ZH} = [ti]{zhi> 


(15) 


where [ti] is the reordered coefficient matrix. Substitution of this equation into 
equation (7) produces the following modified form of the airframe equations: 
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T fr— _ __ " 

[ti] T [mA11][ti] [ti] T [mA12] 


[MA2l][Tl] 


[ti] t [caii][ti] [Tl] X [CA12] 


[CA2l][Tl] 


[KA2l][Tl] 


[CA22] 


[ti] t [kaii][ti] [ti] t [kai2] 


[KA22] 


[MA22] 

flzHlfl 

Jzs} 

/{ziiil 
{zs } J 


fed 

lzs{ 




{ls} 


/ dt = 


To distinguish between the variables of the original and modified finite-element 
models, the variables of equation (16) are renamed: 


f {ZHI J {zh}| 

{zs} I l { zs} / 


= {z} 


And equation (16) is rewritten as 


f Sw a dt = f 

Jo Jo 


6 W a dt = I 6{z} {f} dt = 0 


where 


{f} = [maJ {z} + [ca] { z } + [ka]{z} + {L} 


The matrices [ma] , [ca], [ka], and {l} are defined by comparing 
equation (16). The numbers of variables in {z}, {zh}, and {zs} 
NZH, and NZS . 


equation (19) with 
are denoted by NZ, 


(16) 


(17) 


(18) 


( 19 ) 


The following partitioned form of equation (18) is introduced: 



where 



[MAll][MA12] 

/‘™ } 1 ♦ 

[CAll][CA12] 

[MA2l][MA22] 

\{ZS}J 

[CA2l][CA22] 



+ 


[KAll][KA12] 

[KA2l][KA22] 



( 21 ) 


Formulas for the submatrices in equation (21) are readily discerned by comparison with 
equations (16) and (17). For example. 


m 

[maii] = [ti] [maii][ti] 

(22a) 

{lh} = [tS] t {lh} 

(22b) 


Consistent with equations (7), (8), (21), and (22), Fourier series forms are defined 

for {l}, {lh}, and {ls}: 


{l} = {LOCI + {Lis} sin Qt + {Lie} cos 
{lh} = {LHOC} + {LHls} sin flt + {LHlc} 
{LS} = {LSOC} + {LSIS} sin ftt + {LSIC} 


Qt 

+ ... 

(23a) 

cos 

fit + ... 

(23b) 

cos 

fit + ... 

(23c) 


Definition of an Airframe Trim Solution 

As explained in appendix A, definition of a rotor trim solution provides a basis 
for linearization of the rotor equations. The airframe finite-element model directly 
produces linear equations of motion, so there is no requirement for such a lineariza- 
tion step. However, definition of an airframe trim solution facilitates somewhat the 
derivations shown in this paper. The airframe trim solution is defined to be the 
vector {z} satisfying the following conditions: 

1. The variables in {zh} are constrained to be zero: 


{zh} = {0} 


(24a) 



2. All the variables in {z} are periodic with period 2tt/C2: 


{z (t+2u/£2) } = {Z(t)} 


(24b) 


3. The generalized forces in the vector {fs} given by equation (21) are zero: 


{fs} = {0} 


(24c) 


Equation (24c) expresses that the defined airframe trim solution equilibrates all the 
airframe generalized forces {f} except for the forces {fh} corresponding to the vari- 
ables involved in expressing mechanical connections between the rotor and airframe 
systems. A solution of the finite-element equations meeting the stated conditions 
always exists. 

Notation and partitioning corresponding to the airframe trim solution are indi- 
cated by the following equation: 


{z} trim = {Z0} = 


j{ZHO}| 

kzso}j 


where in view of equation (24a) , 


{ZHO} = {0} 


(25) 


(26) 


The generalized forces corresponding to the airframe trim solution are indicated by 


{P} trim = tF0} ' 



(27) 


where in view of equation (24c ) , 


{fso} = { 0> 


(28) 


The vector partition {zso} which constitutes the unspecified elements of the trim 
solution vector {zo} is determined by the following differential equation obtained 
from equation (21) : 


[MA22]{ZS0} + [CA22]{ZS0} + [kA22]{zSO> + {ls} = {o} 


(29) 
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The generalized forces in the vector {fo} may be interpreted as the external forces 
necessary, in addition to the specified external forces {l}, to maintain the airframe 
variables {z} at the specified trim values {zo}. The vector {fho}, which constitutes 
the unspecified elements of {fo}, is given explicitly by 


{FHO} = [MA12] {ZSO} + [CA12] { ZSO} + [KA12]{ZS0} + {lh} (30) 

which follows from equation (21) . Note that the forces {fho} correspond to the 
variables {zh}, which are the airframe variables directly involved in expressing the 
mechanical connections between the rotor and the airframe systems. Conforming with 
equation (24b) , {ZSO} can be written in Fourier series form as 


{zso} = {zsooc} + {ZSOIS} sin fit + {ZSOIC} cos fit + ... 


(31) 


and it follows that {fho} has the Fourier series form. 


{FHO} = {FHOOC} + {FH01S} sin fit + {FH01C} cos fit + ... 


(32) 


Procedures for computing the airframe trim solution are given in appendices D and E. 


Airframe Equations of Motion Varied From Trim 
The vector of airframe variables {z} is expressed as 


{z} = {zo} + {dz} 


(33) 


where the elements of the vector {dz} represent variations from trim. This equation 
is analogous to equation (Al8a) which was used to linearize the rotor equations. As 
in preceding steps, {dz} is partitioned as 


{DZ} = 


|{dzh}| 

|{dzs}| 


(34) 


From equations (17), (25), and (26), equation (33) may be written as two equations: 


{zh} = {dzh} 


(35a) 


{zs} = {zso} + {dzs} 


(35b) 
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From equation (33), the virtual displacements 6{z} are given by 

6{z} = 6 {dz} (36) 

By using equations (18) , (19) , (33) , and (36) the virtual work equation of an airframe 

finite-element model may be expressed in the form. 


f 6w a dt = f 6{dz} T {{fo} + [ma]{dz} + [ca]{dz} + [ka]{dz}} dt = 0 (37) 

Jo Jo ' 


where 


{fo} = [ma]{zo} + [ca]{zo} + [ka]{zo> + {l} 


(38) 


Equation (37) is equivalent to the differential equations. 


{fo} + [ma]{dz} + [ca]{dz} + [ka]{dz} = {o} 


(39) 


The following partitioned forms of equations (37) and (39) are noted for subsequent 
use in coupling the airframe to the rotor: 



dt 




r 

f{FHO}\ 

[MA11][MA12] 

<< ) + 


1(0}/ 

[MA2l][MA22] 



[CA11][CA12] 

[CA2l][0A22] 



[ka11][ka12] 

[KA2l][KA22] 



dt 


0 


(40) 



[MAll][MA12] 

[MA2l][MA22] 



+ 


[CA11][CA12] 

[CA21][CA22] 



+ 


[KAll][KA12] 

[KA2l][KA22] 




(41) 


Coupling of Rotor and Airframe Equations 

Virtual work for the rotor system and the airframe system is additive. Therefore 
the virtual work equation for the two systems is given by 
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r 



Sw dt 


J (&VJ r + 6w ) dt = 0 
0 


(42) 


Using the partitioned forms which have been established in equations (3) and (40) 
yields for equation (42) 



5W dt = 



r 


AT 


(r 


A 


f{DVBpl 



t { dvh } J 


{ 

\ 

' y*. _ *v 

> l 


^DZH^ 



\{dzs}J 

A 

* 

1 


V 


{qbo} 

{qho} 

{fho} 

.{0}, 




+ 


[MRll][MR12] 
[mR2 1 ] [MR2 2 ] 


[MAll][MA12] 

[MA2l][MA22] 



r A 

f { DVB }") 


\_{dvh}J 

\ 

({dzh^I 


\{dzs}J 
V J 


[CRll][CR12] 

[CR2l][CR22] 

[CA11][CA12] 

[CA2l][CA22] 


r 


a 




C { DVB }! 

\_{dvh}J 

(Jdzs}J 




+ 


[KRll][KRl2] 

[KR2l][KR22] 


[KAll][KA12] 

[KA2l][KA22] 


r 


A 




f {dvb}'1 
\{dvh}J 

f{DZHpl 

\{dzs}J 




J 


dt 


0 


(43) 


The block diagonal form of the three coefficient matrices appearing in equation (43) 
indicates that the rotor equations and the airframe equations are still uncoupled. 
Substituting equations (17) and (35a) into equation (13) yields 


{dvh} = [th]{dzh> 


(44) 


Coupling of the rotor and airframe equations in equation (43) is effected by using 
equation (44) to eliminate the variables {dvh} as independent variables. Formally 
this is accomplished by substituting the transformation. 


L 
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r 




f{DVB}l 


[I] [0] [0] 

r ~\ 

Vdvh}J 

V. J 


[o][th][o] 


{dvb} 

{, ^ 

>- 


< 

{dzh} > 

J{dzh}) 


[0] [I] [0] 


{dzs} 

\{dzs }J 


[o] [0] [l] 


V. J 


( 45 ) 


into equation (43) . This gives the following form for the virtual work equation of 
the coupled system in terms of the independent variables: 



r h T 

{DVB} 


r. 


J <5w dt = I 6<|{dzh}> ( 


{dzs}, 


[MRll] 


[MR12][TH] 
lTr 


[ 0 ] 


[T„] T [MR21] [ ™ ] [” 2 d[TH] [mi2;| 

+ [maii] 


[o] 


[MA21 ] 


[MA22] 


r .. 

{dvb} 




(dzh) / 
{dzs} 


V. 


[CR11] 

[CR12][th] 

[0] 


C • 

{dvb} 

[th] T [cR21] 

[TH] T [CR22][TH] 
+ [CAll] 

[CA12] 

< 

{dzh} 

[0] 

[CA2l] 

[CA22] 


{dzs} 

^ J 

[KRll] 

[KR12][TH] 

[o] 


r 

{dvb} 

[th] T [kR21] 

[TH] T [KR22][TH] 

+ [kaii] 

[KA12] 

< 

{dzh} 

[ °] 

[KA2l] 

[KA22] 


{dzs} 

r 

,1 





{qbo} 


+ < [th] T {QHO> + {FHO} >} dt = 


(46) 


{ 0 } 
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The virtual work equation must be satisfied for any variations of independent vari- 
ables. Therefore, the following differential equation is obtained: 


[mrii] 

[th] T [mR21] 

[o] 


[mR12][th] 

[TH] T [MR22][TH] 
+ [MAll] 

[MA2l] 


[ 0 ] 

[MA12] 

[ma22] 


{dvb} 

< {dzh} > 

{ DZS } 

V > 


+ 


[CRH] 

[TH] T [cR2l] 

Co] 


[CR12][TH] 

[o] 


r "s 

{dvb} 

[TH] T [CR22][TH] 

[CA12] 

« 

{dzh} > 


+ [CAll] 




[cA2l] 

[CA22] 


{dzs} 


+ 


[KRll] [kR12][th] 

[o] 


r 

{dvb} 

[TH] T [KR22][TH] 




[th] T [kR21] 

[KA12] 

< 

{dzh} > 

+ [kaii] 



[o] [KA21] 

[KA22] 


{dzs} 


{qbo} 


<[th] t {qho} + {fho}> 


{ 0 } 


(47) 


Equations (46) and (47) are the virtual work equations and the Lagrange equations of 
motion, respectively, of the coupled rotor-airframe system. 


HARMONIC BALANCE EQUATIONS 

Expansion of the Unknowns in Terms of Rotor Harmonics 

The equations of motion of a coupled rotor-airframe system (eq. (47)) are 
linear differential equations with time-dependent coefficients. The coefficients are 
periodic with period 2TT/Q. The most general periodic solution of these equations 
representing the most general steady-state vibration may be expressed in the follow- 
ing Fourier series form: 


r 

{dvb} 

< {dzh} > 
{dzs} 

V. 


r ~n 

{dvboc} 

<{dzhoc}> + 

{dzsoc} 

v. y 


r 

{dvbis} 

< {dzhis}^ 
{dzsis} 


V. 


y 


r 

{dvbic} 

sin fit + <{DZH1C}> 

{dzsic} 

v y 


COS 


r 

{DVB2S} 

Qt + <{DZH2S}> sin 2Qt 

{DZS2S} 

v. y 


r 

{DVB2C} 

+ <{DZH2C}> cos 2Qt + 
{DZS2C } 

V. y 


r -\ 

{dvbns} 

+ <{dZHNS}> sin NQt 

{dzsns} 

v. y 


r 

{dvbnc} 

/{dzhnc}> COS NQt + 

{dzsnc} 

y 


(48) 
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Following the usual terminology, the first term on the right side in equation (48) 
is referred to as the zeroth harmonic term, the second term as the first harmonic 
sine term, the third term as the first harmonic cosine term, the fourth term as the 
second harmonic sine term, and so forth. When speaking generally, reference is made 
to the Nth harmonic sine and cosine terms. Of course, for N = 0, there is only a 
cosine term, referred to as the constant, or zeroth harmonic, term. 


Reduction of the Virtual Work Equations to Algebraic Equations 


The virtual work equations (eq. (46)) are reduced to a set of algebraic equations 
known as harmonic balance equations by substituting equation (48) into equation (46) 
and carrying out the indicated integration over one rotor revolution. That is, the 
upper limit t is set equal to 2Tr/fi. The requirement that the resulting equation be 
satisfied for independent variations of the unknowns leads to a set of simultaneous 
linear algebraic equations with the following matrix form: 



| {dvboc} j 
, (dzhoc) , 

] {dzsoc} I 

I {dvbis} 1 
/ (dzhis) ) 
{dzsis} 



j {DVB2C } 

< {DZH2C } ) 
I {DZS2C} J 


- 


{qbooc} 

[ [th] t {qhooc} + {fhooc}) 
{ 0 } 

{qbois} 

( [th] t {qhois} + {fhois}^ 
{ 0 } 

{qboic} 

1 [th] t {qhoic} + {fhoic}^ 
( 0 } 

{QB02S } 

[ [th] T {qho2s} + {FH02S} ^ 
( 0 ) 


{QB02C} 

[ [th] t {qho2c} + {fho2c} ' 
{ 0 } 


(49) 

The coefficient matrix in equation (49) is a known constant matrix. Computation of 
this matrix is developed in the next section. The submatrices denoted by hatching 
slanting downward to the right arise from the rotor equations; submatrices denoted by 
hatching slanting upward to the right arise from the airframe equations. In the 
crosshatched areas, contributions from the rotor equations and the airframe equations 
are combined by matrix addition. Blank areas are null. As can be seen, the vector 
on the right side is comprised of contributions from the rotor and airframe trim 
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forces. The reader is reminded that the force vector {FHOOC} is not necessarily 
known because {LOC} in equation (6) is not necessarily known. 


Computation of the Coefficient Matrix of the Harmonic 

Balance Equations 

Airframe contributions .- The airframe contributions to the coefficient matrix of 
equation (49) cluster along the main diagonal of the matrix as indicated by sketch A. 
Except for the zeroth harmonic contributions, these airframe contributions form a 



Sketch A.- Form of coefficient matrix of harmonic balance equations 
with rotor contributions deleted. 
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pattern of uncoupled blocks of four corresponding to the harmonics. Each block of 
four couples the sine and cosine contributions at a given harmonic. The zeroth har- 
monic contributions from the airframe form a single uncoupled block. Expressions for 
the zeroth harmonic block and the Nth harmonic block of four are given by the follow- 
ing equations: 



(50a) 








[ka11][kA12" 

[KA2l][KA22" 

- (Nti) 2 

[MAll][MA12] 

[MA2l][MA22] 


- 

-Nfi 

[CA11][CA12] 

[CA2l][CA22] 

: 

+n£2 

[CAll][CA12] 

[CA21][CA22] 


[kA11][kA12] 

[KA2l][KA22] 

- (Nf2) 2 

[MAll][MA12] 

[MA2l][MA22] 



(50b) 


The need may arise to introduce forces acting on the airframe which are speci- 
fied in terms of impedances. For example, the traditional method of representing 
frequency- independent structural damping leads to forces of this type. Also giving 
rise to such forces are impedance representations of vibration control devices and 
drive system and engine effects. Such forces are introduced through a simple modifi- 
cation of equations (50), as explained in appendix F. 

Rotor contributions .- The rotor contributions to the coefficient matrix of 
equation (49) form an array of blocks as indicated in sketch B. As the sketch shows, 
a four-character indexing system X 2 X 2 X 3 X 4 is used to identify the positions of the 
rotor blocks in the coefficient matrix. In this system, the first character is an 
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Sketch B.- Form of coefficient matrix of harmonic balance equations 
with airframe contributions deleted. 


integer indicating row harmonic, the second character is the letter S or C indicating 
sine or cosine, the third character is an integer indicating column harmonic, and the 
fourth character is the letter S or C indicating sine or cosine. The following 
rule has been derived for computing these rotor blocks: 




[l] 

Co] 

[i] 

[O] 



[dxi x 2 x 3 x 4 ] 


[i] 

[ 0 ] 


[ 0 ] 

[TH]_ 


Co] 

T 

[th] 


[D X 1 X 2 X 3 X 4 11][DX 1 X 2 X 3 X 4 12] 

[D X 1 X 2 X 3 X 4 21][DX 1 X 2 X 3 X 4 22] 


[w 

1 — 1 

0 

1 1 

[[0] 

[TH]_ 


[dx 1 x 2 x 3 x 4 ii] 

[th] T [dX 1 X 2 X 3 X 4 21] 


[ 0 X^ 2 X 3 X 412 ] [th] 

[TH] T [DX 1 X 2 X 3 X 4 22][TH] 


(51) 


where the matrix [dXj_X 2 X 3 X 4 ] is computed according to 


[dx 1 x 2 x 3 x 4 ] = 


(U0K)[KR0C] + ~(ULK) [KRLP] + ^(UHK)[KRHP] 

+ x 3 U ^(uoc) [croc] + |-(ulc) [crlq] + ^-(uhc) [crhq] j 
- (X 3 fi) 2 ^(UOM) [MROC] + i(ULM) [MRLP] + -|-(UHM) [MRHP] ^ 


(52) 


The matrices in equation (52) are the coefficient matrices of the harmonic expansions 
of [mr] , [cr] , and [kr] given in equations (2a) to (2c). Nine integer parameters 
appear in equation (52): UOK, ULK, UHK, UOC, ULC, UHC, UOM, ULM, and UHM. Four 
indices also appear: P, Q, L, and H. These items are computed by an algorithm shown 
in appendix G. The submatrices [0X^X 2 X 3 X 4 ll] , [DX-^X 2 X 3 X 4 12] , [dx^X 2 X 3 X 4 21] , and 

[dX 1 X 2 X 3 X 4 22] appearing in equation (51) are obtained by partitioning the matrix 
[DX 1 X 2 X 3 X 4 ] corresponding to the partitions which are indicated in equation (3) . 

By way of illustration, consider evaluation of the submatrix [d2S2c], In this 
case, X^ = 2, X 2 = S, X 3 = 2, and X 4 = C. Then from the algorithm of appendix G, 

the indices and parameters appearing in equation (52) are 


L = 0 H = 4 P = S Q = C (53a) 

and 


UHK = 1 

UHC = 

1 

UHM = 1 | 

1 


ULK = 0 

ULC = 

0 

ULM = 0 

> (53b) 

UOK = 0 

UOC = 

-1 

UOM = 0 y 



24 



Substituting these quantities into equation (52) yields 


[D2S2C] = — [KR4S] 


<- 


+ 2fi(- [CROC] 


+ — [CR4C] 


)- 


(2fl) 2 ~[MR4S] 


(54) 


Treatment of Airframe Static Forces 

As already noted, to compute the airframe vibrations there is no requirement to 
account for external static forces on the airframe. However, if desired, static 
external forces may be designated to act on the airframe, and such forces are incor- 
porated in the vector {LOC}, as discussed in conjunction with equation (6). Such 
designated static forces are accounted for in the airframe trim solution and are 
eventually reflected in the vector {FHOOC} which is defined by equation (32) and 
appears on the right side of eq uation (49) . As previously noted, there is no require- 
ment for the static forces {LOC} to balance the mean rotor forces represented by 
{qhtr} . As subsequently seen, specification of {LOC} has no effect on the computed 
vibration levels. The only effect is that the vibratory responses are superimposed 
on a static deformation. 

Of course, for a helicopter in steady flight, the static rotor and airframe 
forces must balance. This balance is accounted for by suppressing the mean airframe 
displacements {DZHOC} and {DZSOC}. Thus, the corresponding equations and coefficients 
(rows and columns) are deleted from equation (49) , as indicated by the darkened bands 
in the following equation: 



Ml 

/■ 

[ {DVBIS) 1 


< (dzhis) ) 


{dzsis} . 


I {dvbic} 


< {dzhic} > 


{dzsic} 1 


{ 1 

|{dVB2s} 1 

> = - < 

/ {DZH2S} > 


{DZS2S} 1 


[ {DVB2C} 1 


< {DZH2C } > 


1 {DZS2C} j 


• 


• 


• 

•» 



{qboic} 

[ [tk] t {qhoic} + {fhoic) ^ 
{0} 


{QB02S} 

[ [th] T {qH02S> + {FH02S } ) 
{ 0 } 


{QB02C} 

[ [th] T {qh02C> + {fh02c } ) 
{ 0 ) 


( 55 ) 


25 





Note that the quantities [th]^{QHOOO} and {FHOOC} no longer appear in the equations 
and consequently do not affect the results. If the coupled rotor-airframe system 
admits rigid-body displacements, the removal of the indicated equations and coeffi- 
cients excludes the mean (static) components of such displacements as otherwise 
indeterminate components of the solution. It is very important to recognize that 
oscillatory (N / 0) rigid-body responses are not excluded. 

The assumption that {DZHOC} and {DZSOC} vanish is subject to physical interpre- 
tation. It is equivalent to applying additional static external forces (forces of 
constraint) at the rotor-airframe connection points. These forces may be presumed to 
null the mean displacements {DZHOC}. Then the equations in equation (49) correspond- 
ing to the variables {DZSOC} uncouple with zero right sides; so these variables 
vanish without further addition of external forces. The total forces which must be 
applied at the rotor-airframe interface to suppress the mean interface displacements 
{DZHOC} are designated by {FHTR} and are given by 


{fhtr} = -[th] t {qhooc> 


] 


] 


• • • 


r A 
{dvboc} 

l <°> y 


{ 0 } 

V J 




{dvbis} 

(dzhis) 

{dzsis} 

r ^ 

{dvbic} 

{dzhic}) 

{dzsic} 

r ~\ 

{dVB2s} 

{DZH2S} ) 
{DZS2S } 


{DVB2C } 
{DZH2C } 




{DZS2C } 

v. J 


(56) 
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The row matrix appearing in this equation is comprised of those rows of the coeffi- 
cient matrix of equation (49) corresponding to {DZHOC}. Airframe contributions, 
denoted by hatching slanting upward to the right, have been retained for convenience 
of computation. However, because of the zeros appearing in the column matrix, these 
airframe contributions have no effect on the result. Consistent with the basic 
assumption that the rotor displaces only slightly from the rotor trim solution, one 
can expect that {FHTR} (which should be interpreted as including any forces repre- 
sented by {FHOOC } ) nearly cancels the mean forces (represented by [th] T {qhtr}) speci- 
fied in computing the rotor system trim solution. The cancellation may not be exact 
because in the final solution of the harmonic balance equations, the rotor interface 
variables may undergo oscillatory displacements, whereas these variables are assumed 
to be zero in the trim solution. These additional oscillatory components of displace- 
ment can result in additional increments of mean forces generated by the rotor. 


Alternative Representation of the Rotor System by Impedances 

The rotor contributions to the harmonic balance equations (eq. (49) ) have been 
derived analytically from the linearized equations of motion of the rotor system by 
using the assumed harmonic series given in equation (48) . Alternatively, it has been 
suggested in the literature (see, for example, ref. 8) that the rotor contributions 
to the harmonic balance equations can be obtained by representing the rotor with 
impedances. The impedances would be obtained from numerical integration of the non- 
linear equations of motion of the rotor system. This alternative method of deriving 
harmonic balance equations is discussed in appendix H. 


Definite Forms of the Harmonic Balance Equations 

The harmonic balance equations (eq. (49) ) constitute an infinite set of equations 
with an infinite number of unknowns. To effect a solution in practice, it is neces- 
sary to truncate the general Fourier series for the assumed steady-state solution 
given in equation (48). Also, in equation (48), some variables corresponding to the 
remaining harmonics may be specified to be zero or to have some designated nonzero 
value. The following equation is an example of an assumed definite form of equa- 
tion (48) which might be appropriate for a helicopter with a two-bladed rotor: 


({DVB}) 


r ~\ 

{dvboc} 


<{dzh}) = l 
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{dvbis} 


({dzs]J 


{ 0 } + 


v {0} , 


{ 0 } 
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\ sin fit + L 
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{DVB2S} 

; {o} 

> cos fit + 

{DZH2S } \ 

. (0 » . 
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{DZS2S} 


r -\ 

{DVB2C } 


{DZS2C } 



r 

{DVB3C } 


sin 3^t + < 


{ 0 } 

{ 0 } 


r 

{DVB4S} 


r 

{DVB4C} 


^ cos 3Qt + ^{DZH4S)y sin 4£2t + ^ 
{DZS4S} 


{DZH4C}> cos 4fit 
{DZS4C } 


(57) 
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In this example, all variables corresponding to harmonics greater than four have been 
set equal to zero. The following variables corresponding to harmonics less than four 
have also been set equal to zero: {DZHOC}, {DZSOC}, {dzhIs}, {DZSls}, {dzhIC}, 
{dzsIc}, {DZH3S}, {DZS3S}, {DZH3C}, and {DZS3C}. Specification of zero values for 
the variables (DZHOC) and (DZSOC) may be interpreted as a suppression of static 
rigid- body displacements. Specification of zero values for the remaining variables 
in the list is consistent with the fact that a two-bladed rotor (with perfect balance 
and track) does not exert forces on the airframe at the odd-numbered harmonics. 

In general, dropping all terms beyond a given harmonic involves some approxima- 
tion. Practically speaking, the ordinarily considered vibratory loadings of heli- 
copter airframes decrease very sharply with increasing harmonic number. Therefore, 
useful solutions can be expected with few harmonics included. Suppressing variables 
corresponding to harmonics less than the designated maximum might involve additional 
approximations. However, such specifications are usually based on advance knowledge 
that the specific rotor system does not produce responses with certain harmonics. 

When a definite form is specified for equation (48) , the harmonic balance equa- 
tions (eq. (49)) reduce to a finite number of simultaneous equations with an equal 
number of unknowns. These equations are obtained from equation (49) by setting any 
specified variables to their designated values and then deleting the equations corre- 
sponding to the specified variables. This may be seen by substituting equation (48) 
into equation (46), setting the specified variables to their designated values, and 
restricting the virtual displacements to forms consistent with the designated values 
of the variables. For the example case characterized by equation (57) , equation (49) 
reduces to the following matrix form: 



{dvboc} 

{dvbis} 

{dvbic} 

{dvb2s} 

( {DZH2S} > 

{DZS2S}J 



{DVB3S} 
{DVB3C } 



{qbooc) 

{qbois} 

{qboic} 


I {QB02S} 

UTH] T {0HO2s} + {FH02S}] 

i {0} 


{QB02C} j 

J[th] t {qho2c} + {fho2c}: 
{ 0 } | 


{QB03S } 

{QB03C} 

{QB04S } 

[th] T {QH04s) + {FH04S } 
{ 0 } 

{QB04C} 

[th] t {qho4c} + (FH04c) 
{ 0 } 


( 58 ) 


28 



The finite coefficient matrix is obtained from the infinite coefficient matrix of 
equation (49) by deleting rows and columns corresponding to the specified variables. 

In the example, all specified variables were designated zero. If any of these speci- 
fied variables are designated to have nonzero values, known columns are produced 
which may be transposed to the right. The coefficient matrix is unaffected. Say, for 
example, that some of the variables in the vector {DZH3s} of equation (57) were set 
to nonzero values instead of zero with other specifications of variables in the equa- 
tion unchanged. That is, suppose {0ZH3s} = {d}. Then, the only change to equa- 
tion (58) is an additional column on the right: 



] 

] 

] 


{dvboc} 

{dvbis} 

{dvbic} 



{DVB2C } 

< { DZH2C } \ 

I (DZS2c) I 



Ml 



]M ]HMM ]M ] 



{DVB3S> 
{DVB3C } 



- 


{qbooc} 

{qbois} 

{qboic} 

{QB02S } 

([th] T {qh02S> + {fh02S> ) 
(0} 

{QB02C} 

([th] t {qho2c} + {FH02C}) 
( 0 ) 

{QB03S} 

{QB03C } 

{QB04S} 

( [th] T {qH04S} + {FH04S} ) 
{ 0 } 

{QB04C} 

([th] T {QH04C) +■ {FH04C> ' 
{0} 




{d} 



(59) 


Comments on Instability and Resonance 

The most general solution to differential equation (47) may be expressed as the 
sum of any particular solution and the complementary solution. The steady-state 
solution expressed by equation (48) is a particular solution. The complementary 
solution is defined to be the most general solution of the homogeneous equations 
formed when the vector on the right side of equation (47) is nulled. The comple- 
mentary solution may be stable or unstable. The computational procedures of this 
paper, based on the harmonic balance equations (eq. (49) ) , are not explicitly 
affected by the stability characteristics of the complementary solution. A steady- 
state solution will be generated by the procedures whether or not the complementary 
solution is stable. 
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In vibration analysis, an important consideration is resonance. Resonance is 
indicated when the particular solution increases without bound as some combination of 
system parameters is approached. Resonance conditions are reflected in the harmonic 
balance solution. 


REDUCTION OF HARMONIC BALANCE EQUATIONS USING HARMONIC FORCED 

RESPONSES OF THE AIRFRAME 

The Requirement and the Approach To Reduce the Number of 
Airframe Degrees of Freedom 

The matrices representing airframe contributions to the harmonic balance equa- 
tions (denoted by hatching slanting upward to the right in eq. (49) ) may be large. 
Typical engineering finite-element models of helicopter airframes have several 
thousand degrees of freedom. For practical vibration analysis, the number of air- 
frame degrees of freedom must be reduced. In this paper, this reduction is brought 
about by representing the airframe by forced responses calculated at harmonics of the 
rotor rotational frequency. As subsequently shown, the use of these "harmonic forced 
responses" eliminates from equation (49) all airframe variables except those explic- 
itly involved in expressing the mechanical connections between the rotor and the air- 
frame. The elimination introduces no additional approximations. 

For convenience of discussion, subsequent developments refer to the indefinite 
form of the harmonic balance equations given by equation (49) . From the preceding 
discussion of definite forms, it is straightforward to trace the consequences of 
specifying a definite form of these equations. 


Harmonic Forced Responses 

In this paper, harmonic forced responses for the Nth harmonic (N ^ 0) are defined 
to be a set of column vectors assembled in matrix form as indicated in sketch C: 


[uzhin] 

[uzhon] 

[uzsin] 

[uzson] 

-[uzhon] 

[uzhin] 

-[uzson] 

[uzsin] 


Sketch C.- Matrix of harmonic forced responses for N ^ 0. 
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These forced responses satisfy 


[KAll] [KA12] 


[MAll] [MA12] 

1 

1 

[CAll] [CA12] 


- (Nfi) 2 


' - NS2 

| 


[KA2l] [KA22] 


[MA2l] [MA22] 

1 

1 

[CA2l] [CA22] 



[CAll] 

[CA12] 

1 

I 

[kaii] 

[KA12] 


[maii] 

[MA12] 

+ Nfi 

[CA2l] 

[CA22] 

1 

1 

1 

1 

[KA2l] 

[KA22] 

- (NS2) 2 

[MA2I] 

[MA22] 


[uzhin] 

[uzsin] 

-[uzhon] 

-[uzson] 


[uzhon] 

[uzson] 

[uzhin] 

[uzsin] 


[fuin] 

[fuon] 

[o] 

[0] 

-[fuon] 

[fuin] 

[0] 

[0] 


(60) 

where [uzhon] and [uZSOn] represent out-of-phase responses and [uzhin] and [uZSIn] 

represent in-phase responses. The square submatrices [fuin] and [fuon] may be any 
two matrices for which 


[fuin] [fuon] 
-[fuon] [fuin] 


is nonsingular. Note that this condition is satisfied if either [fuin] or [fuon] is 
nonsingular. It should be recognized that the definition of harmonic forced responses 
herein is a generalization of the ordinary definition of forced response at the fre- 
quency The matrix of harmonic forced responses shown in sketch C is of order 

2(NZ) by 2 (NZH) . Note that only NZH column vectors of this matrix need be computed 
because of repetition in the submatrices comprising the coefficient matrix in 
equation (60) . 

Two commonly used special cases of this general form of the harmonic forced 
responses are the responses to unit imposed harmonic forces and to unit imposed har- 
monic displacements. The first case is obtained by setting [fuin] = [i] and 
[fuon] = [o] in equation (60) . The second case is obtained by setting [uzhin] = [i] 
and [UZHON] = [o]. This can be brought about by a special choice of the forces 
[fuin] and [fuon]. Appendix I discusses this case of responses to unit imposed 
displacements . 

In the case of unit imposed forces, the matrix 

[uzhin] [uzhon] 

-[uzhon] [uzhin] 
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may be identified as the matrix of mobilities of the airframe corresponding to the 
interface degrees of freedom. In the case of unit imposed displacements, the matrix 

[fuin] [fuon] 

-[fuon] [fuin] 


may be identified as the matrix of impedances of the airframe corresponding to the 
interface degrees of freedom. 

Harmonic forced responses for the zeroth harmonic are defined to be a set of 
column vectors assembled in the matrix form indicated in sketch D: 


[UZHIO] 

[uzsio] 


Sketch D.- Matrix of harmonic forced 
responses for zeroth harmonic. 


These zeroth harmonic forced responses satisfy 


[kaii] 

[KA12] 


[UZHIO] 


[fuio] 

[KA2l] 

[KA22] 


[uzsio] 


[0] 


where the submatrix [FUlo] appearing on the right may be any nonsingular matrix. The 
matrix shown in sketch D is of order NZ by NZH. The special case of unit imposed 
forces is obtained by setting [fuio] = [i], and the case of unit imposed displacements 
is obtained by setting [UZHIO] = [i]. 


Computation of Harmonic Forced Responses 

The computation of airframe harmonic forced responses resulting from imposed har- 
monic forces is a standard operation in general-purpose finite-element computer codes 
such as NASTRAN. Thus, the airframe forced responses defined by equation (60) may be 
computed directly. Although it is not necessary for the procedures of this paper, 
there may be some conceptual advantages to utilization of the special forms of the 
harmonic forced responses corresponding to unit imposed forces and unit imposed dis- 
placements. As discussed earlier, responses to unit imposed forces are computed using 
equation (60) after setting [fuin] = [i] and [fuon] = [o]. Direct computation of 
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responses which result from unit imposed displacements ([uzhin] = [i], [uzhon] = [o]) 
is not so convenient with available codes. However, a simple formula for computing 
responses to unit imposed displacements from any solutions of equation (60) is dis- 
cussed in appendix I. 

In design studies, the airframe contributions to the harmonic balance equations 
often must be recomputed many times to study the effects of varying structural members 
and masses or varying impedances of vibration control devices. Appendix E discusses 
a procedure to recompute the harmonic forced responses and the airframe trim solution 
which saves computing effort when only a small number of items are varied. 


Reduction of the Equations 

The reduction of the harmonic balance equations is accomplished by substitutions 
of the following types: 


r ~s 

{dvboc} 


[I] 

< {dzhoc} 

> - 

[UZHIO] 

{dzsoc} 

V J 


[uzsio] 


J {dvboc}! 
UqzhocH 


c > 

{dvbns} 


[I] 

1 

[0] 




{dzhns} 



l 

[uzhin] 1 


[uzhon] 


r 

{dvbns} 

{dzsns} 

- 

_ 

1 

[uzsin] 1 
1 
1 



[uzson] 


{.QZHNS} 

-- - > 

{dvbnc} 


[o] 

1 

Cl] 



{dvbnc} 

{dzhnc} 



1 

-[uzhon] 1 


[uzhin] 


J! QZHNC \ 

{dzsnc} 

V J 



1 

-[uzson] i 


[uzsin] 




(62a) 


(62b) 


These equations may be regarded as transformations of the airframe coordinates by 
which the airframe coordinates are expressed in terms of relatively few new variables 
represented by {QZHOC}, {QZHNS}, and {QZHNC}. Collectively, the transformation equa- 
tions corresponding to the individual harmonics may be expressed as a single trans- 
formation equation which takes the matrix form 
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j { dvboc} 

! {dzhoc} ) 
I {dzsocJj 

|"{dvbis}| 
t {dzhis} ; 
{dzsis} 


| {DVB2C } j 
{DZH2C} 
{dzs2c} 


[I] 


[UZHIO] 

[uzsio] 


[I] 

■ [0] 



[UZHIl] 1 

[UZHOl] 


[uzsil] l 

[uzsoi] 

[0] 

! w 



-[UZHOl] | 

[UZHIl] 


-[uzsoi] | 

[uzsil] 


; [o] 


[UZHI2] ' 

[UZH02] 

[UZSI2] 1 

i 

[UZS02] 

! [i] 


[UZH02] , 

[UZHI2] 

[UZS02] , 

[UZSI2] 


{dvboc} 

{qzhoc} 

{dvbis} 

{qzhis} 

{dvbic} 

{qzhic} 

j {DVB2S } l 
{QZH2S} 
{DVB2C} 
{QZH2C } 


(63) 


Substituting equation (63) in equation (49) leads to equations having the matrix form 



{dvboc} 

{qzhoc} 

{dvbis} 

{qzhis} 

{dvbic} 

{qzhic} 

* {DVB2S } ^ 
{QZH2S} 

{DVB2C } 
{QZH2C } 



(64) 
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There- 


The equations corresponding to the darkened bands are identically satisfied, 
fore, equation (64) reduces to a matrix equation of the form. 



j (dvboc) 1 
KgzHOclj 

j {dvbis)\ 
HqzhisJ j 

J{dvbic}1 
Kqzhic } j 

j {dVB2S} 1 
({QZH2S } j 

j {DVB2C} 1 
){qZH2c } J 


>«- 


{QBOOC} I 

^[th] t {qhooc1 + {fhooc}J 

{qbois} 1 

[th] t (qhois} + {fhois}J 

(oboic} 1 

^[th] t {qhoic} + {fhoic}J 

{QB02s) 1 

[Jth] T {qH02S} + {FH02S}J 

(QB02C) 1 

I [th] T {qH02c} + {fH02c} J 


(65) 


The coefficient matrix in equation (65) is shown as an array of blocks each parti- 
tioned into four submatrices. The difference in shading of the submatrices is for 
convenience in tracking rearrangements of the' coefficient matrix. 


Employing the four-character indexing system defined in a preceding section, the 
typical block in equation (65) is designated [eX-lX 2 X 3 X 4 ] • Symbols for the partitions 
are defined by 



EX 1 X 2 X 3 X 4 


[ex 1 x 2 x 3 x 4 ii] [eX 1 X 2 X 3 X 4 12] 
[eX 1 X 2 X 3 X 4 21] [ex 1 X 2 X 2 X^ 22 "] 


The four submatrices in a typical block are given by 


[EXpX2X 3 X 4 ll] = [0X 1 X 2 X 3 X 4 ll] 


[EX 1 X 2 X 3 X 4 2l] = [th] T [dX 1 X 2 X 3 X 4 21] 


( 66 ) 


(67a) 


(67b) 
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[ex-jX^x^] = [ 0 x^ 2 X 3 X 412 ] [th][uzhix 3 ] 

- a[DX 1 X 2 X 3 Pl2][TH][uZHOX 3 ] 


(67c) 


[EX 1 X 2 X 3 X 4 22] = [TH] T [DX 1 X 2 X 3 X 4 22][TH][UZHIX 3 ] 

- a[TH] T [DX 1 X 2 X 3 P22][TH][uZHOX 3 ] 

+ B[FUIX 3 ] + Y[FUOX 3 ] (67d) 


All matrices appearing on the right sides in equations (67) have been discussed in 
preceding sections. They arise from the rotor contributions to the harmonic balance 
equations, from the coupling equations, and from the airframe harmonic responses. 

The three integer parameters a, 6 / and y and one index P are computed by an 
algorithm shown in appendix J. 

By way of illustration, consider evaluation of the submatrix comprising the 
block [E2S2C] . In this case, X^ = 2 , X 2 = S , X 3 = 2 , and X 4 = C. From the algorithm 
of appendix J, the index and the parameters appearing in equations (67) are given by 

P = S ( 68 a) 


and 


a = -1 8 = 0 Y = 1 

Substituting these quantities into equations (67) yields 


( 68 b) 


[E2S2C11] = [D2S2C11] (69a) 

[E2S2C2l] = [TH] T [D2S2C2l] (69b) 

[E2S2C12] = [D2S2C12][TH][UZHI2] + [D2S2S12] [th] [UZH02] (69c) 

[E2S2C22] = [TH] T [D2S2C22][TH][UZHI2] 

+ [TH] T [D2S2S22][TH][UZH02] + [FU02] (69d) 


As has been discussed, it may be necessary in design studies to vary the airframe 
contributions to the harmonic balance equations (eq. (49)). Of course, this results 
in changes in the coefficient matrix of equation (65) , and a new solution of this 
equation is required. Appendix K discusses a procedure for this reanalysis which can 
save computing effort if the submatrices [EXj_X 2 X 3 X 4 ll] are relatively large compared 
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with the submatrices [EX 1 X 2 X 3 X 422 ] - Some additional comments on solution of equa- 
tion (65) of a purely computational nature relating to order and bandedness of the 
coefficient matrix are given in appendix L. 


RECOVERY OF AIRFRAME RESPONSES 

The total airframe displacements, denoted by the vector {z}, are given by 
equation (33) : 

{ Z } = {zo} + {dz} 

In terms of the partitions defined in equations (17) , (25) , and (34) , this equation 

can be written as 




{dzh}I 

{dzs}, 



{dzh} 


{dzs}I 


(70) 


where {ZHO} has been set equal to zero as in equation (26) . The displacements {zso} 
arc obtained from the airframe trim solution as previously discussed and can be 
written in the Fourier series form shown in equation (31) . The harmonic solutions 
for the displacements {dzh} and {dzs} are obtained by substituting the solution of 
equation (65) into equation (63). These results can be combined to obtain the air- 
frame displacements and the corresponding velocities and accelerations as 


Kzh} 


{dzhoc} 


f {dzhis} 


f {dzhic} 


\{zs } / |{zsooc) + {dzsoc}/ |{zsois} + {dzsis}J 


|{zsoic} + {dzsic}| 


{DZH2C} 


j {DZH2S } 1 j 

< / sin 2 f 2 t + < ? cos 2 !T 2 t + ... 

HZS02S} + {DZS2SM KZS02C} + {dZS2C}| 


1 


(71a) 
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J{z'h} 1 f {dzhis} 1 J {dzhic} 

> = fi< / cos fit - fi\ ) sin fit 

Kzs}J Kzsois} + {dzsis}/ Hzsoic} + {dzsic}) 


{DZH2S} } j {DZH2C} 1 

+ 2fi/ \ cos 2fit - 2fi < S sin 2fit + 

){ZS02S} + {dZS2s}( |{ZS02C} + {dZS2C}| 


(71b) 


J 


J 


{zh} 




{dzhis} 


1 


= - (fiJ^S / sin fit - (fi)' 

Uzs}| /{zsois} + {dzsis}) 


/ 


{dzhic} 




\{zsoic} + {dzsic}/ 


cos fit 


{DZH2S} 1 


{DZH2C} 

- (2fi)^( ) sin 2fit - (2fi)^( ) cos 2fit - ... 

|{zS02s} + {DZS2s}J |{ZS02C } + {dZS2C}| 


(71c) 


If internal member loads are desired, the computed displacements given in equa- 
tion (71a) can be introduced as static displacements imposed on the finite-element 
model of the airframe. For typical finite-element codes, this is a straightforward 
operation. Alternatively, a finite-element code can generate a matrix which, when 
multiplied by {z}, yields desired internal loads. 


The harmonics of the resultant forces acting on the airframe at the interface 
with the rotor are often of interest. These forces are designated by {FZHOC}, 
{FZHls}, {FZHlc}, {FZH2S}, {FZH2C}, and so forth. The forces can be recovered by 
substituting the solution of equation (65) into the right side of 


r ^ 

{fzhoc} 


[fococ] 



r a 

{qzhoc} 

{fzhis} 


[fisis][fisic] 



{QZH1S} 

{fzhic} 


[ficis][ficic] 



{qzhic } 

{FZH2S } 
{FZH2C } 

> - 


[F2S2s][F2S2c] 

[F2C2S][F2C2C] 

< 

{QZH2S} 

{QZH2C} 

V J 



- 


^ 
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where the indicated submatrices, designated [FX-^X 2 X 3 X^], are given by 


[fx 1 x 2 x 3 x^3 = 8[fuix 3 ] + y[fUOX 3 ] 


( 73 ) 


The four-character indexing system and the matrices and parameters appearing on the 
right side of equation (73) are the same as those defined for equation (67d) . 


BASIC COMPUTING SEQUENCE 

With the equations which have been developed, various computing sequences can be 
formed for calculating airframe vibrations. A basic computing sequence is outlined 
here. Two additional sequences are described in appendix M. These incorporate steps 
for reanalysis by the method of appendix E and for representation of the rotor by the 
impedance method of appendix H. 

Figure 4 shows a block diagram indicating the basic sequence of tasks for calcu- 
lating airframe vibrations. Blocks indicating tasks to be performed are arranged in 
three columns to distinguish work done by rotor analysts, work done by airframe 
analysts, and joint work. Supplementary notes follow, keyed by number to the 
individual blocks: 

1. Specify the linearized rotor model : The specification is made by providing 

the coefficient matrices in the Fourier series expansions of [mr] , [cr] , [kr] , and 
{QRO} represented by equations (2) . 

2. jden tify rotor-airframe connection points : The rotor and airframe analysts 

agree on an arrangement of discrete points at which to designate connections. The 
displacements and rotations at these points needed to express connectivity are identi- 
fied. This defines the vector {c} discussed in conjunction with equations (C2) . 

3. Partiti on the rotor equations to isolate interface variables : Among the 

rotor variables {dvr}, a subset {dvh} is identified as the variables which explicitly 
appear in expressions characterizing connections of the rotor model to the airframe 
model (eq. (C2a) ) . The matrices specifying the linear rotor model are partitioned 
accordingly. See discussion of equation (3) and following. 

4. Form roto r contributions to the coupling equations : The matrices [thr] and 

[ger] are formed. See equations (C2a) and (Clla) and related discussion. 

5. Identify harmonics in assumed solution : The steady-state solution given by 

equation (48) is specialized to a definite form agreed on by the rotor and airframe 
analysts. See equation (57) for an example of a definite form. 

6. Specify loads impinging directly on airframe : The distribution of external 

oscillating loads impinging directly on the airframe must be described to enable the 
airframe analysts to define the vector {£} appearing in equation (5) . Static compo- 
nents may be included, but it is not necessary to include them. Loads transmitted 
from the rotor to the airframe through the rotor-airframe connections are not 
included in this specification. 
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ROTOR ANALYSIS 


JOINT ANALYSIS 


AIRFRAME ANALYSIS 



Figure 4.- Block diagram indicating the basic sequence of 
tasks for calculating vibrations. 


7. Form original finite-element model : A finite-element code is employed to 

generate the matrices [ma] , [ca], and [ka], appearing in equation (5). Also generated 
are the load coefficients (lOc), (Lis), (Lie), appearing in equation (6). 

8. Partition the airframe equations to isolate the interface variables ; Among 
the airframe variables {z}, a subset {zh} is identified as the variables which 
explicitly appear in the expressions characterizing connections of the airframe model 
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to the rotor model (eq. (C2b) ) . The matrices characterizing the airframe model are 
partitioned accordingly. See discussion of equation (7) and following. 

9. Form airframe contribution to the coupling equations : The matrix [tha] is 

formed. See equation (C2b) and related discussion. 

10. Form the coupling equations : Compute [Tl] from equation (Clla) . Compute 

[T2] and [T3] from equation (Cllb) . Compute [th] from equation (C14a) . Compute [tc] 
from equation (C14b) . This establishes the coupling equations (9a) and (9b) . 

11. Modify the coupling equations : Identify the submatrices [TCD] and [TCl] 

associated with the matrix [tc]. See equation (10) and related discussion. Identify 
the submatrices [thd] and [THl] associated with the matrix [th]. See equation (12) 

and related discussion. Compute the matrix [th] from equation (14) . The coupling 
equations are thus reduced to a single equation, equation (13) . 

12. Modify the airframe equations : The matrix [Tl] in equation (15) is formed 

from equation (11) , as discussed in the text. Operating with [ti] on the partitioned 
form of the original finite-element equations, as shown by equation (16) , establishes 
the modified form of the airframe equations (eq. (21) ) . Note that this task has been 
to some extent routinized in the NASTRAN code through the feature called multipoint 
constraints. 

13. Compute harmonic forced responses : Equation (60) is solved for each nonzero 

harmonic of interest (see task 5) and the results are assembled as shown in sketch C. 
For any zeroth harmonic responses, equation (61) is solved and the results assembled 
as shown in sketch D. 

14. Compute airframe trim solution : Equation (D3) is solved for {ZHFNS}, {ZSFNS}, 

{ZHFNC} , and {ZSFNC}. The vectors representing trim displacements and forces are com- 
puted from equations (D 6 ) and (D 8 ) . For the zeroth harmonic equations (D9b) to (D9d) 
are used. 

15. Form reduced harmonic balance equations : The matrices [DX-LX 2 X 3 X 4 ] are com- 

puted according to equation (52) and the algorithm of appendix G. Equations (67) and 
the algorithm of appendix J are used to compute the coefficient matrix of equa- 
tion (65). In this process, the only rows and columns generated are those corre- 
sponding to retained harmonics (see task 5) . 

16. Solve reduced harmonic balance equations : Equation (65) is solved. 

17. Recover airframe responses : Equations (63) and (71) are used. If interface 

forces are desired, use equations (72) and (73) . 


DISCUSSION OF THE METHOD 

Two features of the method of this paper seem to have an important bearing on its 
potential utility: (1) use of linearized equations to represent the rotor and (2) use 

of harmonic forced responses to represent the airframe. Comments on the implications 
of these features are offered in the next two sections. 
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Potential Utility of a Linear Rotor-Airframe Model 


With very few exceptions, airframe structural design analysis has always been 
based on linear equations, and this linearity considerably simplifies the analysis. 

It will obviously ease the introduction of vibration analysis into airframe design 
work if rotor models are also linear. A linear rotor model is thus argued to have 
important advantages, but the question remains whether such a model can be adequate. 

The editor of a recent issue of "Vertica" devoted to coupled rotor-fuselage 
dynamics (ref. 20) has commented that even simplified linear models are plausible for 
helicopter vibration analysis: "... evidence available in the literature seems to 

indicate that for aeroelastic response calculations, or the vibration problem of 
coupled rotor/fuselage systems a much simpler linear formulation might be adequate . . . 
The following four reports of correlations of analysis with flight test data also 
bear on the potential of a linear rotor-airframe model: 

1. Reference 8 reports correlations for a tandem rotor helicopter with three- 
bladed rotors. The analysis was based on a rotor-airframe model incorporating linear 
representations of both the rotor and the airframe. The correlations presented are 
reproduced in figure 5. 

2. Reference 10 reports correlations for a compound helicopter with a four-bladed 
hingeless rotor. Plots indicative of the correlations are reproduced in figure 6. 
Nonlinear rotor equations were used directly in the analysis. Therefore, these 
results cannot be cited with complete logic in support of the adequacy of a linear 
rotor-airframe model. The report can, however, be used to evaluate the potential of 

a carefully formed NASTRAN finite-element model, in this case a relatively simple 
model, to represent the airframe response to the major vibration excitations found in 
flight. 


3. Reference 9 reports a correlation for a tandem rotor helicopter with four- 
bladed rotors. The relevant results presented are reproduced in figure 7. The 
analysis was based on a linear rotor-airframe model; however, the linear rotor model 
was very crude. Specifically, only the trim force term {qro} was retained in the 
linearized rotor equations (eq. (A24) ) - the mass, damping, and stiffness terms 
being, in effect, ignored. 

4. Reference 21 reports a correlation performed for a helicopter with a two- 
bladed teetering rotor. In this case, the analysis did not incorporate a model of 
the rotor system. The procedure was to measure the flight vibratory accelerations at 
the rotor hub and then to impose the measured values of acceleration on a NASTRAN 
finite-element model of the airframe. The response of the airframe thus calculated 
was compared with the response measured in flight. Typical results for the major 
responses are shown in figure 8. 

These correlations appear encouraging with regard to the adequacy of linear models 
for coupled rotor-airframe vibration analyses. 

It should be kept in mind that the gross vibratory forces exerted by the rotor 
on the airframe are represented by the vector {QRO} and that this vector would not, 
in general, be computed by linear theory. The vector {QRO} comes from the trim solu- 
tion, and the trim solution, as defined herein, is a solution (by exact or approximate 
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methods) of the underlying nonlinear rotor equations with the constraint that the 
rotor-airframe interface points are fixed. In future structural design work, evalua- 
tions of {QRO} are likely to be strongly influenced by empirical considerations, as 
currently are other airframe design loads such as maneuver loads. The remaining 
terms in the linearized rotor equations (mass, damping, stiffness) are increments of 
the gross rotor forces resulting from displacements of the rotor from the trim solu- 
tion. It is a tenet of design to avoid resonant conditions, and if such conditions 
are avoided, the displacements from trim should be small. Thus, it appears reason- 
able to expect that a rotor model linearized as discussed in appendix A will be 
nearly as good for vibration prediction as the underlying nonlinear model. 
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Figure 5.- Correlation of flight test data with a coupled rotor-airframe 
analysis for 3-per-revolution (3/rev) vibration of a three-bladed 
tandem-rotor helicopter in high-speed level flight. (From ref. 8.) 
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Figure 6.- Comparison of computed coupled rotor-airframe vibrations with measured 
data for compound helicopter with a four-bladed hingeless rotor. Vertical and 
lateral 4-per-revolution (4/rev) vibrations at the forward (STA 131) and aft 
(STA 170) crew stations as a function of airspeed. (From ref. 10.) 


4/rev cockpit vertical vibration 


Vi bratory 
acceleration, 
±g unit 


• o i 
.4 
.3 
.2 
.1 


Flight test 


Analysis 



100 

Airspeed, knots 



Airspeed, knots 


Figure 7.- Correlation of flight test data with an uncoupled rotor-airframe 
analysis for 4-per-revolution (4/rev) vibrations of a four-bladed tandem- 
rotor helicopter. (From ref. 9.) 




Figure 8.- Comparison of flight test data at two flight speeds with 
vibration responses of the airframe calculated using measured hub 
accelerations for 2-per-revolution (2/rev) vertical vibration of 
a two-bladed teetering rotor helicopter. (From ref. 21.) 


Comparison of the Harmonic Forced Response Approach With the Modal 
Approach for Reducing the Number of Degrees of Freedom 
of the Airframe Finite-Element Model 

In an analysis of helicopter vibrations based on a finite-element model of the 
airframe, one must be prepared to reduce the number of degrees of freedom of the 
finite-element model. Two approaches are currently recognized for making this reduc- 
tion and still preserving the essence of the finite-element model: (1) representing 

the airframe by forced responses calculated at a few frequencies corresponding to the 
rotor harmonics of interest, as proposed in this paper, and (2) representing the air- 
frame by superposition of a few of the natural modes of vibration. Whichever approach 
is used, the data (forced responses or modes) needed to represent the airframe with 
relatively few degrees of freedom are calculated by using the finite-element model 
separate from the rotor. 

Modal representations are classic and widely taught. Such representations can 
be used for reducing the number of degrees of freedom when calculating any of the 
linear structural responses which are of interest in practical flight dynamics. This 
includes problems of aeroelastic stability and of transient response as well as the 
present problem of steady-state vibrations. This broad applicability has caused the 
modal representation of the airframe to be the choice of developers of computer simu- 
lations of the helicopter in flight. 

The generality and familiarity of the modal representation is thus acknowledged. 
However, because of the following three advantages, the authors believe that repre- 
senting the airframe system by harmonic forced responses is preferable for vibration 
analyses done to support design of airframe structures: 

1. The analyst can take complete advantage of the fact that the important fre- 
quencies of excitation are usually known ahead of time and are normally few in number. 
Only the forced responses corresponding to excitation frequencies at the important 
harmonics of the rotor rotational frequency need be included in the representation. 
Often it is known in advance that only a single frequency (namely, the so-called blade 
passage frequency which equals the rotor rotational frequency multiplied by the number 
of blades) is important in the loads transmitted from the rotor to the airframe. 

Thus, the analyst has a reliable measure of the adequacy of the reduced-degree-of- 
freedom representation of the airframe. With a modal representation, foreknowledge 

of the important excitation frequencies cannot be comparably utilized. 

2. The computational requirements for reductions based on forced responses should 
be acceptable at all stages of airframe structural design work, including the advanced 
stages when the airframe finite-element model becomes large. This is questionable for 
the modal approach. When computational requirements are the concern, the key differ- 
ence between the two types of reduction is that forced responses are generated by 
solving large sets of simultaneous equations, whereas natural modes are generated by 
solving large eigenvalue problems. Thus, solving for the forced responses can be 
roughly compared with the solutions which are now carried out in sizing the airframe 
for traditional static design loads. Currently the only procedures considered reli- 
able to solve for the airframe natural modes pose much more severe computational 
requirements . 
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3. A reduction based on forced responses lends itself to a relatively quick and 
simple procedure for recomputations when a few selected members of the airframe 
finite-element model are varied in design studies. This advantage also arises from 
the way forced responses are generated, that is, by solution of simultaneous equa- 
tions. As appendix E explains, this allows recomputation by matrix partitioning, a 
well-known technique used in static analyses. This procedure does not result in any 
additional approximations. The only procedures which have been developed for compa- 
rable recomputations of natural modes are very involved and do result in additional 
approximations . 


CONCLUDING REMARKS 

A linear formulation of rotor-airframe coupling intended for airframe structural 
design work has been presented. The airframe is represented by a general finite- 
element model, and the rotor is represented by a set of general linear differential 
equations with periodic coefficients. Coupling of the rotor to the airframe is speci- 
fied through general linear equations. 

Background relating to representation of the airframe is provided by a synopsis 
of finite-element modeling as currently practiced in the U.S. helicopter industry. 

This is supplemented by noting the steps required to incorporate into an airframe 
model special effects which are presented in terms of impedances. Examples of these 
special effects include frequency-independent structural damping, gyroscopic forces 
from rotating engine components, and vibration control devices. 

Background on the rotor representation is discussed to facilitate use of this 
paper by airframe designers. Shown is the process by which linear differential equa- 
tions with periodic coefficients are derived from the usual engineering analysis 
models of rotors. These underlying models are normally expressed in terms of non- 
linear differential equations, which are linearized by assuming small displacements 
from a steady flight trim condition. Also discussed is an alternative representation 
of the rotor suggested in the literature in which the rotor effects are presented in 
terms of impedances at the rotor mounts. 

The coupling equations are derived by assuming that the rotor system and the air- 
frame system are connected at an arbitrary number of discrete points arranged in an 
arbitrary manner. The derivation is general and, in particular, accounts for any 
constraints imposed on the airframe model because of assumptions underlying the rotor 
model. As a supplement, two useful approximate forms of the coupling equations are 
shown. An illustrative example is worked out, along with the derivation of the 
coupling equations, and a logical division of responsibilities (between rotor and 
airframe analysts) for formulation of the coupling equations is suggested. 

The coupling equations are applied to combine the rotor and airframe equations 
into one set of differential equations for the rotor-airframe system. This is accom- 
plished with the displacement method, a customary method for joining components in 
conjunction with the finite-element method of analysis. Factors relating to equilib- 
rium of static components of trim forces are addressed. 

The basis for solving the resulting differential equations to yield the steady- 
state vibrations is the well-known harmonic balance method in which the system steady- 
state vibratory responses are assumed to involve only a limited number of frequencies 
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at integer multiples (harmonics) of the rotor rotational frequency. The paper gives 
an explicit algorithm for generating the harmonic balance equations of a linearized 
rotor-airframe system. 

The number of airframe degrees of freedom is reduced by representing the air- 
frame by harmonic forced responses. These responses are calculated only at the rotor 
harmonics of interest. 

A method is presented for quick recalculation of vibrations when only relatively 
few members of the airframe finite-element model are identified as candidates for 
variation, as often happens in design studies of helicopter vibrations. Constant 
portions of the equations are separated by simple matrix partitioning and solved only 
once; this procedure is well-known in static analysis. Representation of the airframe 
by harmonic forced responses lends itself particularly well to this simple way of 
organizing calculations for rapid calculation of the effects of design changes. 

Explicit relations required to generate the rotor-airframe equations and solve 
for the system vibratory responses are presented in a form suitable for direct com- 
puter implementation. Block diagrams illustrating computating sequences are included 
and discussed. 

Finally, during a discussion of the method, two subjects are addressed: the 

potential utility of linearized rotor-airframe models for practical analysis of 
helicopter vibrations and the merits of representing the airframe by harmonic forced 
responses. It is argued that a properly linearized model should be as good for calcu- 
lating airframe vibrations as the underlying nonlinear model. The procedure of 
representing the airframe by harmonic forced responses is compared with the alterna- 
tive procedure of representing the airframe by superimposing a few of the natural 
modes of vibration. Representation of the airframe by harmonic forced responses is 
recommended for airframe structural design work. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 2, 1982 
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ROTOR SYSTEM EQUATIONS OF MOTION 

This appendix discusses the relationship of linearized rotor equations to funda- 
mental engineering formulations of the rotor system in terms of nonlinear equations. 


Fundamental Equations of Motion 

Engineering analyses of the dynamics of helicopter rotors generally start with 
idealizations which represent the rotor mechanical system in terms of a finite number 
of independent generalized coordinates. These coordinates uniquely determine the 
position of any particle of the system with reference to inertial space. Classical 
techniques of structural mechanics make it possible to derive equations which deter- 
mine the generalized coordinates as functions of time and thereby define the motions 
of the system. These equations are generally nonlinear. The equations allow for 
arbitrary initial conditions and arbitrary external forces. The analyst completes 
the formulation of the fundamental nonlinear equations by specializing the arbitrary 
forces in the equations to represent the forces produced by the interactions of the 
rotor with the airstream. 


Generalized coordinates .- In this paper, the symbols for the generalized coordi- 
nates are Vj , where j = 1, 2, 3, ..., NV. The collection of all generalized 

coordinates is denoted by the vector {v}. As an example, for the rotor blade 
depicted in figure Al, a vector of generalized coordinates is 


r 

K 

6 


M = \ y h } 


z h 

0 


l*J 


(Al) 


The coordinates appearing in equation (Al) are 

x h' Y h' z h fote-and-aft, lateral, and vertical hub displacements with respect to 
inertial axes X^-, Y^, and Zj. 

ip rotation angle of the blade 

0 rigid-body pitch angle 
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3 rigid-body flap angle 

E, single variable multiplying r] (r) , an assigned function of the blade radial 

position r, to describe displacements w(r,t) along the length of the 
blade due to elastic deformation 

As this example illustrates, a generalized coordinate may be defined so that it 
is subject to a direct physical interpretation, or the definition may be abstract. 

In figure Al, the generalized coordinates x h and 4> are in the first category. 

The generalized coordinate E, is an example of an abstract coordinate. Units 
assigned to a generalized coordinate are usually consistent with its definition, such 
as length for the displacement x^. Often a generalized coordinate is regarded as 
nondimensional, such as the rotation angle 1 jj or the coordinate E,. 

Form of equations .- The fundamental equations of motion for the rotor system may 
be presented in the following two equivalent forms: 



(A2) 


2 j 


= 0 


(j = 1, 2, . . . , NV) 


(A3) 


The first form expresses vanishing of a quantity known as virtual work ( 6 w r ) , and the 
second expresses balance of forces. The expressions Qj appearing in these equations 
are usually derived by the methods of Lagrange, and this paper assumes that equa- 
tions (A3) are Lagrange equations. The Lagrange equations are the more familiar of 
the two forms. The virtual work form is included here because it facilitates deriva- 
tion of the coupling of rotor and airframe equations and the introduction of con- 
straints. The quantities 6 vj appearing in equation (A2) are virtual displacements. 

Generalized forces .- The quantity multiplying an increment of a generalized 
coordinate to calculate work is customarily called a generalized force. The quanti- 
ties Oj in equations (A2) and (A3) conform to this definition. It is common and 
quite useful to think of the jth Lagrange equation. 


Qj = 0 (A4) 

as expressing equilibrium of generalized forces. In special cases, an equation such 
as equation (A4) corresponds to equilibrium of forces or torques. For example, when 
a generalized coordinate Vj corresponds to a translation, the expression Qj may 
be interpreted as a force acting in the direction of translation. For the example of 
figure Al, the expression Q 3 would correspond to a force in the fore-and-aft direc- 
tion, and the equation, 

Q 3 = 0 (A5) 

would correspond to equilibrium of external and reaction forces in the fore-and-aft 
direction. Similarly, if the generalized coordinate Vj corresponds to a rotation 
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about some axis, then the expression Qj may be interpreted as a torque about that 
axis. For the example of figure Al, Q 7 would correspond to a torque about the 
shaft axis, and the equation, 

Q 7 = 0 (A6) 

would correspond to equilibrium of applied and reaction torques. 

Relation of general ized forces to generalized coordinates .- Terms in the expres- 
sions Q j arise from two sources : the formulation of the equations of motion of the 

rotor mechanical system and the formulation of the external forces produced by the 
interactions of the rotor with the airstream. For practical purposes, it may be 
assumed that terms arising from description of the rotor mechanical system can be 
expressed as nonlinear functions of the generalized coordinates, their time deriva- 
tives, and time explicitly. The terms arising from description of the aerodynamic 
forces can for many purposes be expressed in the same manner. It is reasonable to 
expect that rotor models used to account for vibrations in airframe design work in 
the future will be relatively simple with the aerodynamic terms so described. This 
renders the generalized forces Qj to be functions of the generalized coordinates, 
their time derivatives, and time explicitly. Any explicit dependence on time is of a 
special character as indicated in subsequent comments on periodicity. Such a rela- 
tionship between the generalized forces and the generalized coordinates is assumed in 
this paper. It is further assumed that time derivatives of the generalized coordi- 
nates higher than the second do not appear in the expressions Q j . Ruling out time 
derivatives higher than the second poses little practical restriction because usual 
mathematical descriptions of rotor systems, including the aerodynamics, are encom- 
passed by the first and second time derivatives of the generalized coordinates. 
Variables with higher time derivatives may be appropriate physically to represent 
certain control system effects in the rotor system equations. However, for linearized 
control system equations, it is straightforward in principle to reduce such a descrip- 
tion to an equivalent description where the highest derivative does not exceed the 
second. Ordinarily, the number of such variables with higher derivatives is quite 
small and, consequently, little computational effort is needed to carry out this 
reduction. 

Prescribed coordinates .- Frequently in applying the fundamental rotor equations, 
it is convenient to designate a specific generalized coordinate v^. to be a pre- 
scribed function of time. Then, in all the expressions Qj the variable Vj. is 
assigned the prescribed form and vj, and v k are obtained by differentiating v k . 

For any coordinate v^ thus prescribed, there can be no virtual displacement. Conse- 
quently, the kth term is deleted from the summation of virtual work (eq. (A2) ) and 
the equation, 

Q k = 0 (A7) 

which is associated with the generalized coordinate v^ is removed from the Lagrange 
equations (eq. (A3)). The expression Qj^, no longer appearing in the Lagrange equa- 
tions, may be interpreted as the generalized force required to constrain v k to the 
prescribed behavior in time. 
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Categorization of generalized coordinates .- It is convenient at this point to 
define four categories of generalized coordinates symbolized by the following parti- 
tions of the vector {v}: 


W 


fr ii 

{vh} 

k - 

W 

{vc} 

( = 1 

{vc} 





(A8) 


The four categories are defined as follows: 

1. Because of the assumption that there is a single rotor, VT is a single 

variable representing the shaft rotation angle. 

2. The vector {vc} represents the variables which are prescribed as functions of 

time to establish steady-state trim. 

3. The vector {vh} represents variables which explicitly appear in expressions 

characterizing mechanical connections of the rotor system to the airframe 
system. 

4. The vector {vb} represents the remaining generalized coordinates of the rotor 

system. 

The vector {vr} represents the elements of {vh} and {vb} collectively. The numbers 
of variables in {vb}, {vh}, {vc}, and {vr} are denoted by NVB, NVH, NVC, and NVR. 

The vector of all generalized forces, denoted by {q}, is partitioned in a corre- 
sponding manner: 


"{QB}" 


m 



l QR f 

{QH} 


L J 

{qc} 

f - s 

{qc} 



V QT, 


(A9) 


The shaft rotation angle VT has been shown separated from the variables in {vh}. 
Recall that {vh} contains the variables explicitly appearing in expressions of con- 
nections between the rotor and airframe systems. This separation arises naturally in 
derivations where displacements of the shaft rotation axis are referred to a coordi- 
nate system fixed in inertial space. In the literature, it is common to reference 
the shaft axis displacements in this manner, but it is also common to reference the 
shaft axis displacements to a coordinate system rotating with the rotor. In the 
latter instance, the separation of VT from the variables in {vh} does not directly 


54 



APPENDIX A 


emerge, but it is straightforward to perform a transformation of variables which 
brings about the separation. 

If there are multiple rotors, the uncoupled equations of motion of the free-body 
rotors may be regarded as a single set of equations. Then the single variables VT 
and QT are replaced by vectors containing the appropriate number of independent rota- 
tion angles and torques. This generalization can be accommodated with only minor 
modifications to the procedures of this paper as long as all rotors are turning at 
the same rotational frequency. Of course, there is a prominent case of multiple 
rotors turning at different frequencies: the helicopter with a single main rotor and 

a relatively small tail rotor to provide antitorque. However, in this case, it is a 
reasonable practice to calculate the airframe vibrations induced by the rotors with 
each rotor treated separately. In multiple-rotor helicopter designs to date where 
the rotors are of comparable size (notably tandem helicopters) , the rotors turn at 
the ■ same frequency. 

Specification of rotor rotational frequency .- The rotational frequency of the 
rotor is introduced by prescribing the rotation angle VT in the form 

VT = Qt + <p (A10) 

where the constant Q is the rotational frequency and the constant (j) is the phase 
angle determining the azimuthal positions of the rotor blades at time equal to zero. 
The equation , 


Q w = o (Ail) 

is removed from the Lagrange equations (eq. (A3)) and the expression may be 

interpreted as the torque required to maintain the rotation at the prescribed 
frequency. 

Comment on periodicity of the nonlinear equations .- Once VT is specified as in 
equation (A10) , the equations of motion of a rotating system have an important prop- 
erty of periodicity. That is, whenever the coordinates Vj are all periodic with 
period 2 7T/Q, then the expressions Qj are all periodic with the same period. 

Stated simply, the generalized forces are periodic whenever the generalized coordi- 
nates are periodic. 


The Trim Condition 

General remarks .- To maintain a specified steady-state flight condition, the 
average forces and moments acting on the rotor must be in equilibrium with the average 
forces and moments acting on the airframe. Such a state of equilibrium is called a 
trim condition. It is important to keep in mind that a trim solution does not gen- 
erally imply constant values of the rotor variables. Some of the rotor variables are 
usually varying (cyclically) . In the analysis of rotor dynamics, the rotor equations 
are commonly linearized by assuming that the rotor motions depart only slightly from 
the trim condition; this assumption has been adopted in this paper. To represent the 
trim condition, the generalized forces {qh} must be maintained with appropriate mean 
values which will be designated {QHTR}. To obtain the mean values requires physical 
interpretation of the coordinates {vh}, knowledge of the rotor-airframe connection. 
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and knowledge of the forces, such as gravitational, inertial, and aerodynamic, acting 
on the airframe. Extensive analysis may be required to make the specification of 
{qhtr}. In this paper, however, it is assumed that the specification has been made. 

Definition of trim solution .- A necessary step in imposing the required mean 
values of {QH} is to compute a trim solution. A trim solution is defined herein to 
be any vector {v} meeting the following conditions: 

1. The variables in {vh} are constrained to be zero: 


{VH> = 0 


(Al2a) 


2. All the variables are periodic in time with period 21T/Q: 




(Al2b) 


3. The mean values of the generalized forces {QH} are 


{QH} 


mean 


{qhtr} 


(A12c) 


Corresponding to any steady flight condition within the operating envelope, there 
should exist an exact solution of the rotor fundamental equations of motion meeting 
the above definition of a trim solution. 

The partitioned vector {v} corresponding to a trim solution is given by 


{v} 


trim 


= {vo} = 


r > 


r A 

{VBO} 


{ VRO / 

{VHO} 

> = l 

l J 

{vco} 

( \ 

{vco} 

fit + 


fit + 


where, in view of equation (A12a) , 


{VHO} = {0> 


(A13) 


(A14) 
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-r- 


The generalized forces corresponding to trim are given by 


E { 2°> 


-< 


c ^ 

{qbo} 


A 

r "i 

{qho} 

> = < 

<QRO> 

{qco} 

/ - \ 

{qco} 



^ T , 


(A15) 


where T is rotor torque and, consistent with the requirement of equation (A12c) , 


{ 2 H0 >mean = ^ HTR > 


(A16) 


Generation of trim solutions .- Various strategies are used to compute trim solu- 
tions (see, for example, refs. 3 and 4) , and an approximate trim solution is commonly 
used in place of the exact trim solution. For instance, trim solutions are sometimes 
obtained by suppressing certain degrees of freedom. For the articulated rotor blade 
of figure Al, for example, one might assume that the blade deformations have a negli- 
gible effect on the trim solution and compute the trim solution with the appropriate 
degree of freedom (£) constrained. 

Variables prescribed for trim .- In arriving at a trim solution, some variables 
(previously denoted collectively by the vector {vc}) are normally prescribed as func- 
tions of time that are not varied in subsequent calculations. For example, in fig- 
ure Al, the blade rigid-body pitch angle 9 might be prescribed to have the form. 


9(t) = 0 Q + 0^ c cos ijj + 0-^ s sin \p 


(A17) 


where 0 Q is the collective pitch angle and 0^ c and 0-^ g are the longitudinal and 
lateral cyclic pitch angles imposed by the flight control system. Equations corre- 
sponding to the variables {vc} are deleted from the Lagrange equations (eq. (A3)), 
the expressions in {qc} being interpreted as the generalized control forces required 
to maintain the specified trim condition. 


Linearization of the Rotor Equations 

The rotor equations of motion are linearized by assuming that vibratory motions 
are small variations (perturbations) from the trim solution. This assumption is 
reflected in the equations 


{VR> = {VRO} + {DVR} (A18a) 

{VC} = {VCO} (A18b) 
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where the elements of the vector {dvr} represent the small variations from trim. 

The generalized coordinates {vc} have been prescribed to establish trim and are not 
varied. As previously discussed, the rotor rotation angle VT is also prescribed 
and is given by equation (A10) . Usual practice is to substitute equation (A18a) for 
the vector of generalized coordinates {vr} into the generalized force expressions Oj 
and to neglect higher order terms in the variables {dvr} occurring in the resulting 
expressions. The following linearized approximation for {qr} is then obtained: 


{qr} = {qro} + [mr]{dvr> + [cr]{dvr> + [kr]{dvr} 


(A19) 


The elements of the square matrices [mr] , [cr], and [kr] are given by 


/3QA 

MR- • = -n- 

13 ' V j/trim 


_ l* Q A 

CR ij \3v i ) . 

' 3 / trim 


KR. . = 


'3Qi' 


i j \ 9v 


3 /trim 


(A20a) 


(A20b) 


(A20c ) 

(i, j = 1, 2, . . . , NVR) 


The subscript trim appearing in equations (A20) indicates that the derivatives are 
evaluated for {v} = {vo}. From the preceding discussion of the periodic character of 
the generalized forces Q j , it can be shown that the matrices characterizing equa- 
tion (A19) are periodic with period 2Tr/fi and therefore may be expressed as Fourier 
series as follows: 


[mr] = [MROC] + [MRls] sin fit + [MRic] cos fit + [MR2S] sin 2fit + . . . (A21a) 
[cr] = [CROC] + [cRls] sin fit + [CRlc] cos fit + [CR2S] sin 2fit + . . . (A21b) 
[kr] = [KROC] + [KRIS] sin fit + [KRlc] cos fit + [KR2s] sin 2fit + . . . (A21c) 
{QRO} = {QROOC} + {QROIS} sin fit + {QROlc} cos fit + (qR 02S> sin 2fit + ... (A21d) 


The coefficient matrices appearing in these equations, such as [MRls] and {QROIS}, 
are constant. 

From equations (A18) and (A10) , the virtual displacements 5 {vr}, 6{vc}, and 6 (VT) 
are given by 
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' a ^nK 
■ 'ft 


6{vr} = 6 {dvr} 

(A22a) 

<5{vc} = { 0 } 

(A22b) 

6(vt) = 0 

(A22c) 


Substitution of equations (A8) , (A9) , (A18) , (A19) , and (A22) into equation (A2) 

yields 


r a t . r 

Jo Jo 


S W r dt = I 6 {dVR> {QR} dt 


= ^ 6{DVR} T |{QRO} + [mr]{dVR} + [cr]{dVR> + [kr]{dVr}| dt = 0 (A23) 


from which it follows that 


{qro} + [mr]{dvr} + [cr]{dvr} + [kr]{dvr} = { 0} 


(A24 ) 


The following partitioned forms of equations (A23) and (A24) are defined for use in 
coupling the rotor to the airframe: 



<5w r dt 



[MRll][MR12] 

[MR2l][MR22] 



[CRll][CRl2] 

f{DVB}l 

[ KR1 1 ] [ KR1 2 ] 

"N 

AdvbA 

[CR2l][CR22] 

\{dvh}/ 

[KR2l][KR22] 

\{dvh}JJ 


(A25) 



[MRll][MR12] 

[MR2l][MR22] 



+ 


[CRll][CR12] 

({dvb}1 

[KRll][KR12] 

[CR2l][CR22] 

|{dvh}J 

[KR2l][KR22] 




(A26) 
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({dvb}\ 

{dvr} = < ) (A27) 

|{dvh}J 


Each of the submatrices appearing in these equations has a Fourier series representa- 
tion similar to equations (A21) . For example. 


[MRll] = [MROCll] 
{qbo} = {qbooc} + 
{qho} = {qhooc} + 


+ [MRlSll] sin fit + [MRlCll] cos fit + . . . 
{QBOls} sin fit + {OBOlc} cos fit + ... 
{OHOIS } sin fit + {QH01C} cos fit + ... 


(A28a) 

(A28b) 

(A28c) 


Equations (A23) and (A24) are the linearized versions of the virtual work equa- 
tion and the Lagrange equations given by equations (A2) and (A3) . As noted in the 
main text, it is assumed that these equations are specified at the outset of the 
computations. The specification is made by providing the coefficient matrices in 
equations (A21) . It is also necessary to specify the variables appearing in {DVH} so 
that the partitions indicated in equations (A25) and (A26) can be made. 
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AIRFRAME SYSTEM EQUATIONS OF MOTION 

This appendix discusses the form of the equations of motion of the airframe 
system. It is assumed that the equations are based on a typical finite-element method 
of analysis. Finite- element methods are used extensively for practical analysis of 
both static and dynamic behavior of flight structures (for commentary, see ref. 22) . 

In particular, the NASTRAN computer code (ref. 1) for finite-element analysis is used 
by most U.S. helicopter airframe manufacturers. 


Finite-Element Models 

Finite-element analysis models are assemblies of structural elements connected 
at discrete points called nodes. Figure Bl and table B1 illustrate a simplified 
finite-element representation of a helicopter airframe. For convenience of illustra- 
tion, the assembly in the figure is drawn as two-dimensional. However, the typical 
assembly would be three-dimensional and may involve many more elements and nodes. 
Structural elements making up these models may include elastic elements in the form 
of rods, shear panels, membrane panels, beam segments, and plate panels, as well as 
concentrated masses and general extended rigid bodies. 

Table Bl, which is to be associated with figure Bl, is a representative tabula- 
tion of information defining the geometry (i.e., topology) of a three-dimensional 
finite-element model. The following comments on the table identify aspects of finite- 
element models which bear on coupling the rotor system to the airframe system: 

1. The numbers in the first column of table Bl identify node points where 
structural elements are joined. 

2. The numbers in the second column designate coordinate systems to which node 
displacements and rotations are referenced. Figure Bl shows two such coordinate 
systems . 

3. The numbers in the third column designate variables associated with each node. 
Up to six variables may be defined. The variables are interpreted as displacements 
and rotations with reference to the indicated node coordinate system. The numbers 1, 
2, and 3 indicate displacements along the X, Y, and Z coordinate axes, respectively, 
and the numbers 4, 5, and 6 indicate the corresponding rotations about these axes. 

A designator missing from this column does not mean that there is a constraint on the 
model. A missing designator means that it is not necessary to define the variable in 
order to formulate the finite-element model. 

4. The designators in the fourth column indicate externally applied forces or 
couples in the direction associated with the designator. 

5. The designators in the fifth column indicate inertial forces or couples 
arising from mass elements incorporated in the model. 

6. The designators in the sixth column indicate suppression of a variable by an 
external constraint posing a restriction on the motions of the model. The example is 
a free-body structure; therefore no designators appear in this column. 
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7. The designators in the last column indicate variables which have specified 
relations to other variables and are therefore not independent. The footnotes to the 
table explain the internal constraints for the example. 

With the exception of the constrained variables designated in the last two columns, 
the variables designated in the third column are independent and may be viewed as the 
generalized coordinates of the airframe finite-element model. 


Zj , 


General Form of the Equations of Motion 

The generalized coordinates of the airframe finite-element model are denoted by 
where j = 1, 2, 3, ..., NZ, or collectively by the vector {z}. As discussed. 


the variables Zj may usually be interpreted as displacements or rotations at nodes. 
As done in appendix A for the rotor system, the airframe system equations of motion 
are first presented in general virtual work and differential equation forms which 
are equivalent : 


r^-n 

Jo Jo jTi 


6Zj 


F j dt 


(Bl) 


F . = 0 
3 


(j = 1, 2, 


NZ) 


(B2) 


The symbols Fj, which are denoted collectively by the vector {f}, represent the 
generalized force expressions for the airframe. Consistent with the interpretation 
of the generalized coordinates, the generalized forces may usually be thought of as 
net forces and couples at nodes. For the prevalent finite-element methods of analysis, 
the generalized force expressions have the form. 


{f} = [ma]{z[ + [ca]{z} + [ka]{z> + {£} 


(B3) 


Insertion of equation (B3) into the matrix forms of equations (Bl) and (B2) gives the 
following forms for the virtual work equation and the differential equations of the 
airframe finite-element model: 


J <Sw a dt = J <5{z} T |[ma]|z} + [ca]{z| + [ka]{z} + {l}| dt = 0 


(B4) 


[ma] { z i + [ca]{z} + [ka]{z> + {L> = {0} 


(B5) 


The square matrices [ma], [ca], and [ka] are constant, real, and symmetric and are 
either positive definite or positive semidefinite. They are the mass, damping, and 
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stiffness matrices. These matrices may be quite large, since the number of general- 
ized coordinates may range up to several thousand; but the matrices are generally 
sparsely populated. The vector {l} represents the external forces acting on the air- 
frame. These forces in general may be arbitrary functions of time. 

The following partitioned form of equation (B4) is defined for use in coupling 
the rotor to the airframe: 



Svtf dt 

cl 




[MAll][MA12] 

[MA2l][MA22] 



[CAll][CA12] 

[CA2l][CA22] 



where 


+ 


[KAll][KA12] 

[KA2l][KA22] 



dt 


0 



(B6) 


(B7a) 


(B7b) 


The vector {zh} represents variables which explicitly appear in expressions charac- 
terizing mechanical connections of the rotor system to the airframe system. The 
vector {zs} represents the remaining variables of the airframe finite-element model. 
The numbers of elements in {zh} and {zs} are denoted by NZH and NZS . 


Comment on Application of Linear Constraints to the Airframe Finite-Element 


Model Implied by Assumptions Made in Forming 


the Rotor Model 


As discussed in the text and in appendix C, coupling 
airframe model may imply additional linear constraints on 
model. Such constraints take the form 


the rotor model to the 
the airframe finite-element 


[c]{z> = {0} 


(B8) 


Note that equation (9b) is a special case of this equation. Application of such con- 
straints is a standard procedure in finite-element modeling of structures. For the 
NASTRAN finite-element code, in particular, inclusion of linear constraints is 
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systematic as long as the analyst can identify a nonsingular square submatrix of [c] 
of order equal to the number of rows of [c] (see eq. (9b) and related discussion) . 

If the rows of [c] are linearly independent as in equation (9b) , at least one such 
submatrix always exists. 

With regard to the practical problem of identifying an appropriate nonsingular 
square submatrix, the following conditions are expected when coupling rotor models to 
airframe models : 

1. The number of equations of constraint and the number of variables involved 

are not large. 

2. The coefficients of the equations normally can be presented in algebraic form. 

3. The engineer is normally working with a complete physical interpretation of 

the equations. 

Under these circumstances, the required identification of a nonsingular square sub- 
matrix should not prove difficult. It is noted, however, that on strictly numerical 
grounds, that is, when the matrix [5] is presented in numerical form only, systematic 
identification of such a nonsingular square submatrix may be difficult. 
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TABLE Bl.- REPRESENTATIVE INFORMATION DEFINING GEOMETRY OF A 
THREE-DIMENSIONAL FINITE-ELEMENT MODEL 


Node 

point 

Coordinate 

system 

Node point 
degrees of freedom 
(a) 

External 

forces 

(a) 

Inertial 

forces 

(a) 

External 

constraints 

(a) 

Internal 

constraints 

(a) 

1 

1 

123 

_ 

123 

_ 

_ 

2 

1 

123 

3 

- 

- 

- 

3 

1 

123 

- 

123 

- 

- 

4 

1 

123 

- 

- 

- 

- 

5 

1 

123 

3 

123 

- 

- 

6 

1 

123 

- 

123 

- 

- 

7 

1 

123 

- 

123 

- 

- 

8 

1 

123 

3 

123 

- 

- 

9 

1 

123 

- 

123 

- 

- 

10 

1 

123 

- 

123 

- 

- 

11 

1 

123 

- 

123 

- 

- 

12 

1 

123 

- 

123 

- 

- 

13 

1 

123 

- 

123 

- 

- 

14 

1 

123456 


- 

- 

b 123456 

15 

1 

123456 

123456 

123456 

- 

- 

16 

1 

123456 

- 

- 

- 

b 123456 

17 

1 

123456 

- 

123456 

- 

- 

18 

1 

123456 

- 

- 

- 

b 123456 

19 

1 

123 

- 

123 

- 

c 123 

20 

1 

123456 

- 

123456 

- 

- 

21 

1 

123 

- 

123 

- 

c 123 

22 

1 

123 

- 

123 

- 

- 

23 

1 

123 

- 

123 

- 

- 

24 

1 

123 

3 

- 

- 

- 

25 

1 

123 

- 

- 

- 

- 

26 

1 

123 

- 

- 

- 


27 

1 

123 

3 

123 

- 


28 

1 

123 

- 

123 

- 

- 

29 

1 

123 

- 

123 

- 

- 

30 

1 

123 

3 

- 

- 

- 

31 

1 

123 

- 

- 

- 

- 

32 

1 

123 

- 

- 

- 

- 

33 

1 

123 

3 

123 

- 

- 

34 

1 

123 

- 

123 

- 

- 

35 

1 

123 

- 

123 

- 

- 

36 

2 

123 

- 

123 

- 

- 

37 

2 

123456 

- 

123456 

- 

d l 

38 

2 

123456 

2 

123456 

- 

d l 

39 

2 

123456 

2 

123456 

- 

d l 


a The numbers 1, 2, and 3 designate displacements along X, Y, and Z coordinate 
axes; 4, 5, and 6 designate rotations about these axes. 

^Related to degrees of freedom of transmission. 

c Related to degrees of freedom of engine. 

d Axial deformation of vertical tail beam suppressed. 
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COUPLING EQUATIONS 
Connection Points 

It is assumed that the rotor system and the airframe system are joined at a 
finite number of discrete points. For illustration, figure Cl shows a simplified 
coupled system. In this two-dimensional example, there are two connection points 
denoted by the open circles. Connection points may, but do not necessarily, coincide 
with nodes of the airframe finite-element model, which are denoted in the drawing by 
solid circles. Associated with each of the connection points, variables representing 
up to three displacements and three rotations are defined as described for node points 
in appendix B in the section on finite-element models. The displacements and rota- 
tions, for all connection points, are denoted by Cj, where j = 1, 2, ..., NI, or 

collectively by the vector {c}. 


Implied Constraints 

Before coupling equations are established, it is important to recognize that 
because of different assumptions in the modeling of the rotor system and the modeling 
of the airframe system, connecting the two models together may imply constraints on 
either of the systems. Three different situations can occur: 

1. A constraint is implied on the airframe model, as would occur in the example 

of figure Cl if the transmission were modeled as rigid. 

2. No constraints are implied, as would occur in the example if the transmission 

were modeled as flexible. 

3. A constraint is implied on the rotor model. 

Situations (1) and (2) are addressed by this appendix. Situation (3) is unlikely and 
is not addressed. However, the required modifications to the computational procedures 
to account for this situation are straightforward. 


Basic Relations 

According to the definition of {vh} and {zh} as vectors of generalized coordi- 
nates, it is always possible to establish relations connecting {vh} to {c} and {zh} 
to {c}. These relations are represented by the following two equations: 


{f r } 

= {c} 

(Cla) 


= {c} 

(Clb) 


where the elements of the vector {f r } are functions of the variables VHj 
(j = 1, 2, ..., NVH) , and the elements of the vector { f a J are functions of the 
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variables ZHj (j = 1, 2, . .., NZH) . The functions in {f r } and {f a } are generally 

nonlinear. 

The rotor displacements have been assumed to vary only slightly from a trim 
solution, as indicated by equation (A18a) . Because the variables in {vho} for a trim 
solution are by definition equal to zero, the variables in {vh} are approximated by 

the small variation variables {dVh}. Consistent with usual assumptions in finite- 

element analysis, the variables in {zh} are assumed to be small. It follows that 
equations (Cl) can be approximated by linear forms indicated by 


[thr]{dvh} = {c} 

(C2a) 

[tha]{zh} = {c} 

(C2b) 


The matrices [thr] and [tha] are constant matrices. For practical interface arrange- 
ments, the matrix [thr] has linearly independent columns and the matrix [tha] has 
linearly independent rows. This implies that the number of rows in [thr] is equal to 
or greater than the number of columns and that the number of columns in [tha] is 
equal to or greater than the number of rows . 


Figure C2 illustrates specification of the basic relations. In this two- 
dimensional example, which shows the rotor and airframe models separated for pictorial 
clarity, {vh} is a three-element vector where VH^ and VH 2 are hub displacements and 
VH^ is the hub pitching rotation. In this case {c} is a four-element vector repre- 
senting displacements at the two connection points. With the rotor shaft and trans- 
mission assumed to be rigid, equation (Cla) becomes 



sin 

vh 3 

+ 

%2 

(1 " 

COS 

VH 3 f 


r 

c l 

VH 2 + £ 2 

sin 

VH 3 

+ 

£ 1 

(1 - 

cos 

vh 3 ) 

\ = 

c 2 

VH 1 - £ 1 

sin 

VH 3 

- 

£ 2 

(1 - 

cos 

V h 3 ) 

? \ 

c 3 

^VH 2 - * 2 

sin 

vh 3 

+ 

£ 1 

(1 - 

COS 

vh 3 )^ 


^ c 4 - 


Assumption of small values for the variables on the left of equation (C3) leads to a 
linear form corresponding to equation (C2a) : 


— 


— 


r 

1 

0 

-&1 


r 


C 1 





DVH-^ 



0 

1 

%2 



V 1 

C 2 




< 

DVH 2 

> = < 


1 

0 

-*1 




c 3 





DVH 3 



0 

1 

_o 

£ 2_ 




^ C 4 > 


(C4) 
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To illustrate another way for specifying the basic relations, equation (Clb) is 
determined by linear interpolation: 


a 



ZI^ 



a 





c 2 


b 

£4 


zh 3 



c 3 







(C 5 ) 


Equations (C 5 ) are already linear and therefore no linearization step is required. 
The equations take the following matrix form corresponding to equation (C 2 b) : 


a 


A 



—— 


r ■> 


f 

Cl 

£3 

0 

1 " h 

0 

0 

0 


ZE ± 


C 1 








z 5 2 



0 

cL 

£3 

0 

1 " £3 

0 

0 

J 

zh 3 

> -1 

c 2 



b 


, b 


\ 

zh 4 

f \ 


0 

0 

£4 

0 

1 - — - 
£4 

0 




c 3 








ZH5 



0 

0 

0 

b 

0 

1 


zh 6 


c. 




£4 


£4 

v > 

< 4 7 


(C6) 


Derivation of Coupling Equations 
The basic relations (eqs. (C 2 ) ) are first written in the form 



{ 0 } 


(C 7 ) 
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For example, equations (C4) and (C6) written in this form are 


1 0 -A-jl 

0 1 l 2 

1 0 -SL 


0 


1 




'o' 


0 

V. J 


(C8) 

Note that the number of individual equations represented by equation (C7) is NI where, 
as previously defined, NI is the number of variables in the vector {c}. The use of 
only elementary row operations reduces equation (C7) to. 


[Tl] | -[T2] 
[0] ' -[T3] 




(C9) 


where [Tl] is square and nonsingular and the rows of both [t 2] and [t 3] are linearly 
independent. For the example, this is accomplished by interchanging rows 3 and 4 and 
then subtracting row 1 from row 4. These elementary row operations yield 



(CIO) 
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which is of the form of equation (C9) . Such a reduction can be carried out syste- 
matically by one of the various forms of Gaussian elimination using only elementary 
row operations. Formally any such Gaussian elimination procedure applied to equa- 
tion (C7) can be expressed as the following two equations: 


[ger][thr] 


[Tlf 

[ 0 ] 


(Clla) 


[ger][tha] 


[T2] 

[T3] 


(Cllb) 


where [ger] is a square nonsingular matrix reflecting the elementary row operations. 
For the example. 


[ger] = 


1 

0 

0 

0 


1 

0 

0 

0 


1 

0 

0 

0 

0 

1 

0 

0 


0 

1 

0 

0 


0 

1 

0 

0 

0 

0 

1 

0 


0 

0 

0 

1 


0 

0 

0 

1 

-1 

0 

0 

1 


0 

0 

1 

0 


-1 

0 

1 

0 


(C12) 


In cases where coupling the rotor and airframe equations implies constraints on the 
airframe equations, a logical division of responsibilities would appear to be the 
following: 

1. The rotor analyst provides [thr] and [ger] 

2. The airframe analyst provides [tha] 

with both working to a commonly accepted definition of the connection points and the 
connection variables {c}. 

Equation (C9) is now written as two equations: 


{dvh} = [th]{zh} 


(Cl 3a) 


[tc]{zh} = {o} 


(C13b) 



APPENDIX C 


where 


[th] 

= [Tl] 1 [T2] 

(C14a) 

[tc] 

= [T3] 

(C14b) 


The numbers of individual equations represented by equations (C13a) and (C13b) are NVH 
and (NI - NVH), respectively, where, as previously defined, NVH is the number of rotor 
variables which explicitly appear in the relations connecting the rotor to the air- 
frame. The rows of both [th] and [tc] are linearly independent. 

In the special case where the number of variables in {dvh} equals the number of 
variables in {c}, the matrix [thr] appearing in equation (C7) is square and non- 
singular. In this case, no airframe constraints as represented by equation (C13b) are 
defined, and the equation defining the matrix [th] which appears in equation (C13a) 
follows directly from the basic relations as 


[th] = [thr] 1 [tha] 


(C15) 


As has been shown, equations (C13) are equivalent to the linearized basic rela- 
tions given by equations (C2) . Equation (C13a) expresses the rotor interface vari- 
ables {dvh} explicitly in terms of the airframe interface variables {zh}. Equa- 
tion (C13b) involves only the airframe interface variables and can be recognized as 
constraints on the airframe finite-element model inherent in the formulation of the 
rotor model. Equations (C13) are adopted in this paper as the general form of the 
coupling equations . 


Approximate Formulation of the Coupling Equations 

As stated, equations (C13) follow exactly from the basic relations given by 
equations (C2) . As a practical matter, engineers sometimes formulate coupling equa- 
tions which are more or less inconsistent with the basic relations. Such formulations 
may ignore the basic relations altogether as when fictitious massless structural ele- 
ments are incorporated to establish coupling, or the basic relations may be taken into 
account in an approximate way. Two such approximate procedures are noted here. 

In the first procedure, certain of the constraints represented by equation (C13b) 
are ignored. This is permissible if it can be assessed that the airframe is con- 
structed so that it behaves approximately as specified in some or all of the constraint 
equations. For the example of figure C2, the constraint equation (C13b) becomes the 
single equation. 
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5-3 


0 



0 




ZH 2 


- 3 

ZH 4 

zh 5 

ZHf- 
o J 


(C16) 


Assume for simplicity that = & 4 = £ and that a = b = SL/2. The equation then 
reduces to 


ZH 


1 


= ZH 5 


(C17) 


That is, the horizontal displacements at the two finite-element nodes indicated in 
figure C2 are equal. Since these nodes are spanned by axial members, it might very 
well be judged that the pertinent stiffnesses are large, so that equation (C17) is 
approximately satisfied. This procedure may often be justified. However, it is 
judicious to derive equation (C13b) and evaluate separately the implications of 
retaining or discarding each one of the individual constraint equations. 

In the second procedure equation (C7) is approximately satisfied on a least 
squares basis. Explicit equations of constraint corresponding to equation (C13b) do 
not emerge. The coupling relation corresponding to equation (C13a) is obtained by 
minimizing the quadratic form, 


){dVHH , -i T r f{ DVt 

) [[thr] j -[tha]J |Jthr] i -[tha]J < _ 


{dvh}[ 

{zh} 


with respect to the rotor variables {DVH} which yields 
{dvh} = [Jthr] t [thr]J [thr] t [tha]{zh> 


(C18) 


Equation (C18) is in the form of equation (C13a) where 


[th] = |^[thr] t [thr]J [thr] t [tha] 


(C19) 


This least squares procedure may be useful in preliminary design where the details of 
the mechanical connection between the rotor and the airframe are not yet defined. 
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CALCULATION OF AIRFRAME TRIM SOLUTION USING 
AIRFRAME HARMONIC FORCED RESPONSES 

The differential equations defining the nonzero portions ( {zso} , {FHO}) of the 
airframe trim displacement vector and the trim force vector are given by equations (29) 
and (30) . Writing equation (21) in harmonic balance form while taking into account 
equations (29) , (30) , (26) , and (28) leads to equations of the following form to be 

solved for the harmonic components (N ^ 0) of the trim solution: 


[KAll] [KA12] 
[KA2l] [KA22] 


(nS 3) 2 


[maii] 

[MA2l] 


+ NSJ 


[CAll] 

[CA2l] 


[CA12] 

[CA22] 


[MA12] 

[MA22] 


[kaii] 

[KA21] 


- 


[CAll] 

[CA2l] 


[CA12] 

[CA22] 


[KA12] 

[KA22] 


(NJ2) 


2 [MAll] 
[MA21] 


[MA12] 

[MA22] 


r *\ 



r to -i 


Hfhons} - {lhns}] 

KzsonsH 


\ - {lsns} j 

/ N 

> = < 

[ 

/ ) 


Hfhonc} - {lhnc }1 

W ZSONC } J 


[ - {lsnc} j 


t J V J 


(Dl) 


For the purposes of this appendix, the solution for the displacement vector is now 
written in the form. 


r "\ 

{0} 


r n 

{ZHFNS} 


[uzhin] 

[uzhon] 

{ZSONS} 


{ZSFNS> 

k + 

[uzsin] 

[uzson] 

_ 

{0} 

\ 

{ ZHFNC } 

? + 

-[uzhon] 

[uzhin] 

{ZSONC} 
V. J 


{ZSFNC} 

V J 


-[uzson] 

[uzsin] 



(D2) 


where the vectors {ws} and {wc} are as yet undetermined and the first vector on the 
right side is obtained from solution of 


[KAll] [KA12] 
[KA2l] [KA22] 


(N«) 


2 [MAll] 
[MA2l] 


[MA12] 

[MA22] 


- NS}' 


[CAll] 

[CA2l] 


[CA12] 

[CA22] 


{ZHFNS} 

{ZSFNS} 


" A 

{lhns} 

{lsns} 


-------- 


+ Nn 


[CAll] 

[CA21] 


[CA12] 

[CA22] 


[KAll] [KA12] 
[KA2l] [KA22] 


(Nf2) 


2 [mah] 

[MA2l] 


[MA12] 

[MA22] 


{ZHFNC} 

{ZSFNC} 

V 


{lhnc} 

(lsnc) 


(D3) 
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Note that the harmonic forced responses defined in conjunction with equation (60) and 
sketch C appear in equation (D2) . The vectors {LHNS}, {lsns}, {lhnc}, and {LSNC} 
appearing on the right of equation (D3) are associated with any oscillatory forces 
impinging on the airframe from external sources, and these vectors are assumed to be 
known. The vectors {ws} and {wc} are chosen to satisfy 


J {ZHFNS} 1 

1 {zhfnc } J 


[uzhin] 
- [uzhon] 


[uzhon] 

[uzhin] 




which yields 



[uzhin] 

-[uzhon] 


[uzhon] 

[uzhin] 


■ i u 


{zhfns} 1 
zhfnc } J 


(D4) 


(D5) 


Substituting equation (D5) into equation (D2) then gives the nonzero portions of the 
airframe trim displacement vector as 


J {zSONS} 1 

Uzsonc } J 


j {ZSFNS}! 


|JZSFNC}J 


[uzsin] 

-[uzson] 


[uzson] 

[uzsin] 


[uzhin] 

-[uzhon] 


[uzhon] 

[uzhin] 


-1 


f{ zhfns}! 
! {zhfnc}/ 


(D6) 


The corresponding airframe trim force vector is given by 


{fhons} 


[fuin] 

[fuon] 

{0} 

\ _ 

[0] 

[0] 

{fhonc} 

/ - 

-[fuon] 

[fuin] 

. {0} ^ 


[0] 

[°] _ 


f(ws} 

i(wc) 


VrV . 


(D7) 


which yields 


J{fhons} 

Ufhonc} 


[fuin] [fuon] 
-[fuon] [fuin] 


[uzhin] [uzhon] 
-[uzhon] [uzhin] 


-1 


J {zhfns} 

|^{ ZHFNC }^ 


(D8) 
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As pointed out in the main text, it is not necessary to account for the static 
external forces acting on the airframe. However, if desired, static external forces 
may be designated to act on the airframe and can be accounted for by allowing zeroth 
harmonic components in the trim solution. The needed zeroth harmonic equations corre- 
sponding to equations (Dl) , (D3) , (D6) , and (D8) are 


[kaii] 

[KA2l] 


[kaii] 

[KA2l] 


{ZSOOC} = {ZSFOC} - [uZSIO] [UZHI0] _1 {ZHF0C} (D9c) 

{FHOOC} = -[FUI0][UZHI0] _1 {ZHF0C} (D9d) 


[KA12] 

[KA22] 


f {0} 

1 {zsooc} 


{fhooc} - {lhoc}| 
- {lsoc} 


(D9a) 


[KA12] 

[KA22] 


J{ZHF0C} 
( {ZSFOC} 



(D9b) 


The vector {FHOOC} appears in the harmonic balance equations (see eqs . (49) and (65)). 

However, note that if the procedures leading to equation (55) are used, the equations 
containing {FHOOC} are deleted; therefore {FHOOC} has no effect on calculation of the 
vibratory responses. The only effect of static forces on the airframe is to super- 
impose a static displacement on the vibratory responses (see eq. (33)). 
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RECOMPUTATION OF AIRFRAME HARMONIC FORCED RESPONSES 
AND AIRFRAME TRIM SOLUTION BY MATRIX PARTITIONING 

In design studies, the need often arises to recompute the airframe contributions 
to the harmonic balance equations to study the effects of varying structural members 
and masses or varying the impedances of vibration control devices. This appendix 
formulates a computational procedure based on the use of simple matrix partitioning 
to recompute both the harmonic forced responses of the airframe and the airframe trim 
solution. When relatively few masses and structural members are candidates for varia- 
tion considerable savings in computing effort can be realized from this procedure. 
Reference 23 is a recent report on the same subject. 

For purposes of this appendix, the airframe forced response equations for the 
Nth harmonic (N ^ 0) given in equation (60) are written in the more general parti- 
tioned form. 



The partitioning in equation (El) reflects categorization of the airframe variables 
as follows: 

1. The first and fourth block rows correspond respectively to the sine and 

cosine components of the interface variables {dzh}. The number of equations 
in each of these block rows is NZH. 

2. The second and fifth block rows correspond respectively to the sine and 

cosine components of variables which define the deflections of the elements 
which are varied. The number of equations in each of these block rows is 
NZSV. 

3. The third and sixth block rows correspond respectively to the sine and cosine 

components of the remaining variables. The number of equations in each of 
these block rows is NZSF. 

It is recognized that to categorize variables in this manner and to cast equation (60) 
into the form of equation (El), in general, requires renumbering of variables and 
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reordering of rows and columns. For simplicity of discussion, no rearranging is 
assumed. When rearranging is necessary, the appropriate procedures should be 
apparent . 

As an illustration of this categorization of airframe variables, consider the 
simplified finite-element model of a helicopter airframe given previously in fig- 
ure Bl. Assume that the interface with the rotor (rotor not shown) is the single 
node numbered 15 which corresponds to the rotor hub. Assume further that masses and 
structural members to be varied are distributed throughout the airframe as depicted 
in figure El. Then the degrees of freedom associated with node point 15 would be in 
the first category, the degrees of freedom associated with node points 1, 11, 14, 18, 
19, 20, and 30 to 38 would be in the second category, and the degrees of freedom 
associated with the remaining node points would be in the third category. 

Equation (El) can be reduced to the following two equations: 



[All] 

[A12] 

-[bii] 

-[B12] 


[A2l] 

[A22] 

— [B2l] 

-[B22] 


[bii] 

[B12] 

[All] 

[A12] 


[B2l] 

[B22] 

[A2l] 

[A22] 


[FI] 


[A13] 

-[B13] 

[F2] 


[A23] 

-[B23] 

[Gl] 


[B13] 

[A13] 

[G2] 

_ _l 


[B23] 

[A23] 


[A13] 

[A23] 

[B13] 

[B23] 


-[B13] 

-[B23] 

[A13] 

[A23] 


[A33] 

-[B33] 

“1| 

[A3l] 

[A32] 

- [b31] 

-[B32] 

[B33] 

[A33] 

! 

[B3l] 

[B32] 

[A3l] 

[A32] 


[A33] 

-[B33]~ 

-1 

[F3] 

[B33] 

[A33] 


[G3] 


[XI] 

[X2] 

[Yl] 

[Y2] 


(E2a) 


[xsf 


[A33] 

-[B33] 

-1 

[A3l] 

[A32] 

-[B31] 

-[B32] 


[XI fj 


[A33] 

-[B33] 

-1 

[F3] 

[Y3] 


[B33] 

[A33] 


[B3l] 

[B32] 

[A3l] 

[A32] 


[X2] 

+ 

[B33] 

[A33] 


[G3] 










[Yl] 















[Y2] 







(E2b) 


Equation (E2a) is solved for [xi], [x2], [Yl], and [Y2], and the result is substituted 

into equation (E2b) to obtain [X3] and [Y3]. The coefficient matrix in equation (E2a) 
is relatively small (order equal to twice the total number of degrees of freedom in 
categories 1 and 2) , but contains the inverse of a large matrix (order equal to twice 
the number of degrees of freedom in category 3) . However, the matrix which must be 
inverted is not affected by varied masses and structural members, and thus the inverse 
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need be computed only once. In this connection, it should be emphasized that while 
this matrix is formally indicated as an inverse matrix in equations (E2) , in practice 
matrix decomposition techniques would be employed to compute matrix products involving 
an inverse as in these equations. These decomposition techniques require much less 
computational effort than does formal inversion. The effect of varied masses and 
structural members is contained solely in the submatrices [A22] and [B22], which 
appear only in the left term of the coefficient matrix in equation (E2a) . 

To compute the harmonic forced responses of the airframe using equations (E2) , 
the following identifications are made: 


[XI] 


[uzhin] 

[uzhon] 

[X2] 


[UZSIN2] 

[UZS0N2] 

[X3] 


[UZSIN3] 

[UZS0N3] 

[Yl] 


-[uzhon] 

[uzhin] 

[Y2] 


-[UZS0N2] 

[UZSIN2] 

[Y3] 


-[UZS0N3] 

[UZSIN3] 

[FI] 


[fuin] [fuon] 


[F2] 


[o] 

[o] 


[F3] 


[o] 

[0] 


[Gl] 


-[fuon] [fuin] 


[G2] 


[0] 

[0] 


[G3] 


[0] 

[0] 


[uzsin] = 

[UZSIN2] 

[UZSIN3] 




(E3a) 


(E3b) 


(E4a) 


[uzson] 


[UZSON2] 

[UZSON3] 


(E4b) 
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As a basis for computing the solution for the airframe trim displacement vector, 
equation (D2) must be formed. Equation (E3a) provides the ingredients of the coeffi- 
cient matrix of the second term on the right of equation (D2) . To compute the first 
column on the right of equation (D2) , the following identifications are made in 
equations (E2) : 


[XI] 


r ' 

{ZHFNS} 

[X2] 


{ZSFNS 2> 

[X3] 


{ZSFNS 3} 

[Yl] 

- \ 

{ ZHFNC } 

[Y2] 


{ZSFNC2} 

[Y3] 


{ZSFNC3} 
V. J 


(E5a) 


[fi] 


{lhns} 

[F2] 


{LSNS 2 } 

[F3] 


{LSNS3> 

[gi] 

_ 

~ \ 

{lhnc} 

[G2] 


{LSNC2> 

[G3] 


(LSNC3) 

V. J 


(E5b) 


where 


{ZSFNS} 


J {ZSFNS2}] 
1 {ZSFNS3} J 


V, 


J 


(E6a) 


{ZSFNC} 


{ ZSFNC2} 1 
{ ZSFNC 3 } J 


(E6b) 


{lsns} = 


J {LSNS2} l 
HLSNS3} J 


(E6c) 
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{LSNC3 } J 
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The previous two computations can be combined by identifying 


[XI] 


[uzhin] [uzhon] {zhfns} 

1 — 1 
CN 

X 
1 1 


[UZSIN2] [UZS0N2] {ZSFNS2} 

[X3] 


[UZSIN3] [UZS0N3] {ZSFNS3} 

[yi] 


-[uzhon] [uzhin] {zhfnc} 

[Y2] 


-[UZS0N2] [UZSIN2] {ZSFNC2} 

[Y3] 


-[UZSON3] [UZSIN3] {ZSFNC3} 

[FI] 


[fuin] [fuon] -{lhns} 


[F2] 


[o] [0] -{LSNS2} 


[F3] 


[0] [0] -{LSNS3} 


[gi] 


-[fuon] [fuin] -{lhnc} 


[G2] 


1 1 

0 

1 i 

i — i 

0 

1 1 

K 

S 


[G3] 


[O] [0] -{LSNC3} 



For the zeroth harmonic, the equations corresponding to equations (El), 
and (E7) are 


[All] [A12] [A13] 


[yi] 


[Gl] 

[A2l] [A22] [A23] 


[Y2] 

= 

[G2] 

[A3l] [A32] [A33] 


[Y3] 


' r—| 

0 
U) 

1 1 

1 


(E6d) 


(E7a) 


(E7b) 


(E2) , 


(E8) 
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[All] 

[A12] 


[A13] 

[A33] 1 [[A3l] [A32]] 


[yi] 


[gi] 


[A13] 

[A33] 1 [G3] 


[A2l] 

[A22] 


[A23] 



[Y2] 


[G2] 


[A23] 



(E9a) 


[Y3] = -[A33] X [[A3l] 



[yi] 

[Y2] 


+ [A33] 1 [G3] 


(E9b) 


[YI] 


[UZHIO] {ZHFOC} 

[Y2] 

= 

[UZSI02] {ZSFOC2} 

[Y3] 


[UZSI03] {ZSFOC 3} 

[Gl] 


[fuio] -{lhoc} 

[G2] 


[0] - {LSOC2 } 

[G3] 


[0] -{LSOC3} 


where 


[uzsio] 


[UZSI02] 

[UZSI03] 


(ElOa) 


(ElOb) 


(Ella) 


{ZSFOC } = 


J {ZSFOC2}! 
1 {ZSFOC3 }J 


(Ellb) 


{lsoc} = 


j {LSOC2}! 
1 {lsoc 3 } ( 


V. 


(Ellc) 


Coordinate 
system 1 



Figure El.- Masses and structural members of airframe finite-element model which are to be varied. 
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AIRFRAME FORCES WHICH ARE PRESENTED IN TERMS OF IMPEDANCES 


The need may arise to introduce forces acting on the airframe which are specified 
in terms of impedances. For example, the traditional method of representing frequency- 
independent structural damping leads to forces of this type . Also giving rise to such 
forces are impedance representations of vibration control devices, effects of drive 
system rotation, and effects of engine rotation. Forces represented by impedances are 
introduced directly into the harmonic balance equations given in equation (49) . Let 
these additional forces be designated by the vector {fdz}. The Fourier series expan- 
sion of this vector, in partitioned form, may be written as 


{FDZ} = 


{fdzh}1 J {fdzhoc} 1 J {fdzhis} 1 

{fdzs}j 'Wfdzsoc}J Y{fdzsis}j 


sin fit + 


J {fdzhic} 1 
1 {fdzsic } J 


cos Qt 


I {FDZH2S} 1 

1 (fdzs2s ) J 


sin 2Qt + 


J {FDZH2C } 1 
1 {FDZS2C } J 


cos 2flt + . . . 


(FI) 


where 


f { fdzhoc}"} 

{{ FDZSOC } J 

f {FDZHIS}) 
\_{FDZS1S}/ 

({FDZHIC}') 
{{ FDZS1C } j 

{{FDZH2S}') 
ljFDZS2S } J 

{ C{ FDZH2C }) 
{{ FDZS2C } j 




HIO 


— - 


- 

HIl 


HOI 

HOI 


HIl 

_ _ 


_ __ 





(’{dzhoc}'! 

\Jdzsoc}J 
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Ijozsis}/ 

f{DZHlc}'l 

\{dzsic}J 

{{DZH2S}') 
\{DZS2S } J 

f {DZH 2 c}) 
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(F2) 


f{FDZHNS}') 

\j>dzsns}J 

f{FDZHNC}'l 

(Jfdzsnc}/ 
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and where [hio] is the constant (zeroth harmonic) impedance and [hin] and [hon] are 
the in-phase and out-of-phase Nth harmonic (N ^ 0) impedances. Equation (F2) is in 
the nature of a collection of impedance relationships expressing force amplitudes as 
linear functions of displacement amplitudes. The typical (Nth harmonic) impedance 
equation in this collection is given by 


r 

J {fdzhns }1 
Ufdzsns } J 

\ _ 


HIN 

HON 

_ _j 


J 

r > 

j {dzhns}1 

) {dzsns }J 

J {FDZHNC }1 

) {fdzsnc}J 
V ^ J J 

f - 

- 

HON 

HIN 


\ 

J {dzhnc } 1 

| {dzsnc}j 
V ^ J J 


The coefficient matrix in equation (F3) is known as an impedance matrix and is deter- 
mined by the linearized characteristics (mechanical, electrical, or otherwise) of a 
particular device. It is a square matrix, constant in time and generally a function 
of the rotational frequency The block diagonal form of equation (F2) , which 

uncouples the harmonics, is appropriate for devices represented by linear equations. 

The effect of the resulting additional airframe forces on equation (49) is 
reflected in the following modifications to equations (50) : 







[KAll][KA12] 

[KA2l][KA22] 
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[hio] 
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[KA11][KA12] 

[KA2l][KA22] 
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[MA2l][MA22] 


-Nfi 

[CA11 ] [CA1 2] 

[CA2l][CA22] 






+Nfi 
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- (NO) 2 
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[MA2l][MA22] 







" 




[hin] 



[hon] 
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' 







-[hon] 



[hin] 
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For the special case where the additional forces represent frequency-independent 
structural damping, the matrix [hIn] is null and the matrix [hOn] is independent of 
the rotational frequency A. 
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ALGORITHM FOR COMPUTING THE PARAMETERS 
AND INDICES IN EQUATION (52) 

The formula given in equation (52) for the matrix [dx 3 X 2 X 3 X 4 ] contains the nine 
integer parameters - UOK, ULK, UHK, UOC, ULC, UHC, UOM, ULM, and UHM - and the four 
indices - P, Q, L, and H. The algorithm for computing these items is given in this 
appendix. Some of the conditions for computing the parameters are not mutually 
exclusive. In such cases, the last condition always takes precedence. 

Computation of indices 

L = |x ] _ - X 3 | 

H = X x + X 3 


P = c 

Q = s 

(if 

x 2 = x 4 ) 

P = s 

to 

II 

n 

(if 

x 2 ^ x 4 ) 


Computation of parameters for stiffness terms 


UHK = 



r 

-1 



(If 

X 2 = s 

and 

X 4 = S) 

- ) 

i 



(If 

x 2 = c 

and 

X 4 = C) 

\ 

i 





(If 

*2 * X 4> 


L 0 



(If 

x 1 = 0 

and 

x 3 = 0) 


r 

1 





(If 

x 2 = x 4 ) 


-1 

(If 

x 2 = S 

and 

X 4 = c 

and 

x 3 > x x ) 


1 

(If 

x 2 = c 

and 

x 4 = s 

and 

x 3 > x x ) 

- 

1 

(If 

X 2 = s 

and 

X 4 = c 

and 

x 3 < x x ) 


-1 

(If 

x 2 = C 

and 

X 4 = s 

and 

x 3 < x 4 ) 


0 





(If 

X 1 = x 3) 


0 





(If 

X 1 = 0) 
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UOK = 


Computation of parameters for damping terms 


UHC = < SN 

lo 


ULC = 


UOC = 


(If 

X^ = X 3 and 

x 2 = X 4 ) 

(If 

X x ? X 3 or 

x 2 ^ X 4 ) 

terms 




(If 

x 4 = S) 


(If 

X 4 = C) 


(If 

X 2 * X 4 ) 


(If 

X 2 = X 4 ) 


(If 

x 3 = 0 ) 


(If 

X 2 ? x 4> 

(If 

X 2 = X 4 and 

X 3 > Xl ) 

(If 

X 2 - X 4 and 

*3 < x l) 


(If 

X 1 = X 3 ) 


(If 

X 

H 

II 

O 


(If 

X 3 = 0) 

(If 

X-^ = X 3 and 

X 2 * X 4> 

(If 

X 1 ^ X 3 or 

X 

11 

tN 

X 


(If 

X 1 = 0) 


(If 

X 3 = 0) 


Computation of parameters for mass terms 


UHM = UHK 
ULM = ULK 
UOM = UOK 


(If X, jt 0) 


UHM = 0 ULM = 0 


(If X 3 = 0) 
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REPRESENTATION OF THE ROTOR SYSTEM BY IMPEDANCES 


In this paper, the rotor contributions to the harmonic balance equations (see 
sketch B and eqs. (51) and (52)) have been derived analytically on the basis of 
linearization of the rotor equations and harmonic expansion of the unknown rotor dis- 
placements (see eqs. (1) and (49)). It has been suggested in the literature (see, 
for example, ref. 8) that it would be feasible to formulate the rotor contributions in 
the form of impedances computed directly by numerically integrating the fundamental, 
generally nonlinear, equations of motion of the rotor system. This appendix explains 
this approach. 

Imagine the rotor system trimmed in a steady flight condition. Assume that the 
rotor is idealized and represented by fundamental nonlinear equations as discussed in 
appendix A. Then the rotor variables involved in the rotor-airframe interface are 
denoted by {vh}, and the corresponding generalized forces required to maintain the 
trim condition are denoted by {qh}. If the rotor displacements are assumed to vary 
only slightly from a trim solution as discussed in appendix A, then the variables {vh} 
may be represented by the small motion variables {dvh} which have the following 
Fourier series form: 


{dvh} = {DVHOC} + {DVHls} sin ftt + {DVHlc} cos Qt + 


(HI) 


The total forces represented by (QH) and the trim solution forces represented by {QHO} 
may be expressed in Fourier series forms as 


{QH} = {QHOC} + {QH1S} sin Qt + {QHIC} cos fit + ... (H2) 

{QHO} = {QHOOC} + {QH01S} sin fit + {QHOlc} cos fit + ... (H3) 

Under the assumption of small displacements, the coefficients {QHOC}, {QHls}, and so 
forth in equation (H2) may be expressed by a linear equation of the form. 
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{qhois} 
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{dvhis} 
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L J 

{qhoic} 


[kuicoc][kuicis][kuicic][kuic2s][kuic2c] . . . 


{ DVHlc} 

{QH2S} 

f = 1 

{QH02S} 

) + 

[KU2S0C][KU2S1S][KU2S1C][KU2S2S][KU2S2C] . . . 

< 

{DVH2S} 

{QH2C} 


{QH02C} 


[KU2C0C][KU2C1S][KU2C1C][KU2C2S][KU2C2C] . . . 


{DVH2C} 

V " ) 


V ' ^ 


• • * 
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Equation (H4) expresses the force-displacement relationship for the rotor system in 
terms of the interface degrees of freedom. The vector of trim solution forces given 
by the first term on the right side of equation (H4) and the submatrices [kUX^X 2 X 3 X 4 ] 
which comprise the coefficient matrix in equation (H4) can be computed as described in 
reference 8 by numerically integrating the nonlinear equations of motion of the rotor 
system. It should be recognized that these quantities must in principle be derived 
anew whenever the flight condition changes because the validity of equation (H4) is 
predicated on linear behavior of the rotor system in the neighborhood of a trim 
solution. 

Equation (H4) has been set forth to determine the generalized forces {qh} repre- 
senting external forces required to maintain specified generalized displacements {dvh}. 
Alternatively, the formulation may be viewed as equations, in harmonic balance form, 
which determine the displacements {dvh} corresponding to specified generalized 
forces {qh}. The harmonic components of the specified generalized forces appear in 
the left column. In the special case for which the left column is nulled, as in the 
following equation, the equation represents a free-body rotor system with no forces 
at the hub: 


[kuococ][kuocis][kuocic][kuoc2s][kuoc2c] . . . 


r A 

{dvhoc} 


r 

{qhooc} 

[kuisoc][kuisis][kuisic][kuis2s][kuis2c] . . . 


{dvhis} 


{qhois} 

[kuicoc][kuicis][kuicic]Ckuic2s][kuic2c] . . . 


{dvhic} 


{qhoic} 

[kU2S0c][kU2S1s]CkU2S1c][kU2S2s][kU2S2c] . . . 

< 

{DVH2s} 

> --< 

{QH02S} 

[kU200c][kU2C1s][kU2C1c][kU2C2s]CkU2C2c] . . . 


{DVH2C} 


{QH02C} 

• 


V. • > 


V • J 


Equation (H5) can be specified at the outset of computations in lieu of equa- 
tions (2) to represent the rotor system. In this event, the harmonic balance equa- 
tions (eq. (49) ) have the form 
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(dzhoc) 


r >1 

[th] t {qhooc> + (fhooc) 

{dzsoc} 


{0} 

(dzhis) 


[th] t {qhois} + {fhois} 

(dzsis) 


{0} 

(dzhic) 


[th] t {qhoic) + {fhoic } 

(dzsic) 


(0) 

1 {DZH2S) 


[th] T {QH02s} + {FH02s} f 


{DZS2S} {0} 


{DZH2C } 
{DZS2C} 


[th] T {QH02C} + {FH02C} 
{ 0 } 


J 


(H6) 


Equation (51) for computation of the rotor contributions to the coefficient matrix is 
replaced by 



= [th] t [kux 1 x 2 x 3 x 4 ][th] 


(H7) 


Note that in this approach the variables {DVB} do not appear. 

Further, when the rotor system is represented by impedances in this manner, 
the transformation equation (63) takes the form 
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f {dzhoc}'! 
Hdzsoc} J 


[UZHIO] 

[uzsio] 


r { dzhis} 
[Jdzsis} 
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f{DZHlCpl 

\{dzsic}J 
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f{ DZH2C} 
[_{DZS2C } 


-[UZH02] [UZHI2] 
-[UZS02] [UZSI2] 


r a 

{qzhoc} 

{qzhis} 

{qzhic} 

{QZH2S} 

{QZH2C} 






J 


(H8) 


and the structure of the reduced harmonic balance equation (65) becomes 


[][][][][]• • 


r *\ 

{qzhoc} 


”■ "\ 
[th] t {qhooc} + {fhooc} 

[][][][][]• • • 


{qzhis} 


[th] t {qhois} + {fhois} 

[][][][][]• • - 
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[][][][][]• • • 
[][][][][]• • • 
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{QZH2C } 
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[th] T {QH02S} + {FH02S} 

> 

[th] T {QH02C} + {FH02C } 

• • • • • • 

• • • « • • 


<- > 




(H9) 
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The submatrices of the coefficient matrix in equation (H9) are denoted by [gx^X 2 x 3 x 4] 
and are given by 

[gx.,x 2 x 3 x 4 ] = [th] t [kux 1 x 2 x 3 x 4 ][th][uzhix 3 ] 

- a[TH] T [KUX 1 X 2 X 3 P][TH][uZHOX 3 ] 

+ B[FUIX 3 ] + y[fUOX 3 ] (H10) 


which is a specialization of equation (67d) . As discussed in connection with equa- 
tion (67d) , the parameters a, (3, and y and the index P appearing in equation (H10) 
are computed by rules given in appendix J. 

The total airframe displacements , and the corresponding velocities and accelera- 
tions, are recovered with equations (71) . When impedances are used to represent the 
rotor as discussed in this appendix, the harmonic solutions for the airframe dis- 
placements {dzh} and {dzs} are obtained by substituting the solution of equation (H9) 
into equation (H8) . These results, and the airframe trim solution, are then substi- 
tuted into equations (71) for the final recovery of the displacements, velocities, 
and accelerations. If needed, the resultant forces acting on the airframe at the 
interface with the rotor are obtained by substituting the solution of equation (H9) 
into equation (72) . 
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CALCULATION OF HARMONIC FORCED RESPONSES FOR 
UNIT IMPOSED HARMONIC DISPLACEMENTS 

Finite-element computer codes for computing harmonic forced responses are gen- 
erally predicated on imposed harmonic forces rather than on imposed harmonic displace- 
ments. However, if needed, harmonic forced responses associated with unit imposed 
displacements are readily obtained from the harmonic forced responses calculated for 
imposed forces. Using the harmonic forced responses resulting from imposed forces for 
the Nth harmonic (N ^ 0) , one can define 


[f] 

[g ] 1 


[uzhin] 

[uzhon] 

-[g] 

[f] 


-[uzhon] 

[uzhin] 


(ID 


Postmultiplying equation (60) by the left matrix in equation (11) results in response 
and force matrices of the forms shown in figure II. The harmonic forced responses 
corresponding to unit imposed displacements are given in figure II (a) and the corre- 
sponding set of forces required to maintain the unit displacements are given in 
figure II (b) . 

The corresponding matrices for the zeroth harmonic are obtained by defining 


[f] = [UZHIO] 


-1 


( 12 ) 


and postmultiplying equation (61) by [f]. This yields the responses and forces shown 
in figure 12. 
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[I] 

[o] 

[uZSIN][f] 

[uZSIN][g] 

- [lJZSON][g] 

+ [uZSON][f] 

[0] 

[l] 

-[uZSON][f] 

-[uZSON][g] 

- [uZSIN][g] 

+ [uZSIN][f] 


(a) Responses. 


[fuin] [ f ] 

[FUIN] [g] 

- [FUON][g] 

+ [fuon] [ f ] 

[o] 

[o] 

-[FU0N][f] 

-[FU0N][g] 

- [FUIN][g] 

+ [fuin] [ f ] 

Co] 

[°] 


(b) Forces. 

Figure II.- Harmonic forced responses corresponding to 
unit imposed displacements calculated from responses 
corresponding to imposed forces. 


[l] 


[fuio] [f ] 

[uzsio][f] 


[°] 


(a) Responses. (b) Forces. 

Figure 12.- Zeroth harmonic responses corresponding to 
unit imposed displacements calculated from responses 
corresponding to imposed forces . 
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ALGORITHM FOR COMPUTING THE PARAMETERS AND INDEX IN EQUATIONS (67) 

The formulas given in equations (67) for the submatrices into which matrix 
[ex^X 2 X 3 X 4 ] is partitioned contain the three integer parameters a, 3r and y 
the index P. The algorithm for computing these items is as follows: 


a = 0 (If X 3 = 0) 


a = 

1 

P = c 





(If 

x 4 = S) 

a = 

-1 

P = s 


(If 

X 3 

/ 0 

and 

x 4 = C) 

6 = 

1 

o 

II 

>- 


(If X 

1 = 

X 3 

and 

X 2 = X 4> 

3 = 

0 

rH 

II 

(If 

X-^ = X 3 and 

X 2 

= S 

and 

x 4 = C) 

3 = 

0 

Y = -1 

(If 

— X 3 and 

X 2 

= C 

and 

X 4 = S) 

3 = 

0 

o 

II 

>~ 





(If 

X x / X 3 ) 


and 
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SOLUTION OF REDUCED HARMONIC BALANCE EQUATIONS BY MATRIX PARTITIONING 

As pointed out in the main text and in appendix E, the need often arises in 
design studies to recompute the airframe contributions to the harmonic balance equa- 
tions in order to assess the effects of varying structural members and masses or 
varying impedances of vibration control devices. This means that the reduced harmonic 
balance equations given in equation (65) must be solved each time the airframe struc- 
ture is changed. This appendix formulates a computational procedure based on the use 
of simple matrix partitioning to solve the reduced harmonic balance equations. In 
recomputation, considerable savings in computing effort can be realized from this pro- 
cedure when the rotor representation involves many degrees of freedom relative to the 
number of degrees of freedom required to describe the interface with the airframe. 

To this end, equation (65) is first rearranged as follows: 


li 

II 

if 


'/s'/// S' 

is / ' U Xls 


1 Mil II 
I il 11 li 
Will I! 

mmm 

111 MB 


t 

rinr 

e ][ it 

Hill 
ill ill 


II li 
Bill 

III li 

] 
1 


{dvboc} 

(dvbis) 

{DVBIC} 

{DVB2s} 

(DVB2C) 


(qzhoc) 

(qzhis) 

(qzhic) 

{QZH2S} 

(QZH2C) 


{qbooc} 
{qbois} 
(qboic) 
(QB02S) 
{ QB02C } 


< ) 

[th] t {qhooc) + (fhooc) 

[th] t {qhois} + (fhois) 

[th] t {qhoic) + {fhoic} 

[th] T {qh02s) + (FH02S) 

[Th] T {QH02c) + (FH02C) 


(Kl) 
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Let equation (Kl) be represented by the equation 


M ; [b] 

_tc] ; [ D ] 

Equation (K2) can be reduced to the following two equations: 
[[D] - [c][ A r 1 [B]]{x2} = {P2} - [C][A] _1 {P1} 

{xi> = [a] -1 {{pi} - [b]{x 2 }} 



{PI} 


{P2} 


(K2) 


(K3a) 

(K3b) 


Equation (K3a) is solved for {x2} and the result substituted into equation (K3b) to 
obtain {xi}. The coefficient matrix in equation (K3a) is relatively small (maximum 
order equal to (NT) (NZH) by (NT) (NZH) ) but contains the inverse of a relatively large 
matrix (maximum order equal to (NT) (NVB) by (NT) (NVB) ) . However, the matrix which 
must be inverted does not contain the effects of airframe changes and thus the inverse 
need be computed only once. As noted in appendix E, matrix decomposition techniques 
would be used rather than formal inversion. The effects of airframe changes are con- 
tained solely in the submatrices [b] and [d]. The computations needed to evaluate the 
matrix products involving [b] can be reduced in situations requiring recomputation by 
writing [b] in the form 


[b] = [Bl][El] + [B2][E2] 


(K4) 


where 


[bi][ei] = 


[DOCOC12][DOC1S12][DOC1C12][DOC2S12][DOC2C12] . . . 
[D1S0C12][D1S1S12][D1S1C12][D1S2S12][D1S2C12] . . . 
[d1C0C12][d1C1S12][d1C1C12][d1C2S12][d1C2C12] . . . 
[D2S0C12][D2S1S12][D2S1C12][D2S2S12][D2S2C12] . . . 
[d2C0C12][d2C1S12][d2C1C12][d2C2S12][d2C2C12] . . . 


[th][uzhio] 

[th][uzhii] 

[th][uzhii] 

[TH][UZHI2] 

[TH][UZHI2] 


(K5a) 


101 


APPENDIX K 


and 


[B2][E2] 


[□OCOC12][DOC1C12][DOC1S12][DOC2C12][DOC2S12] . . . 
[D1S0C12][D1S1C12][D1S1S12][D1S2C12][D1S2S12] . . . 
[D1C0C12][D1C1C12][D1C1S12][D1C2C12][D1C2S12] . . . 
[D2S0C12][D2S1C12][D2S1S12][D2S2C12][D2S2S12] . . . 
[D2C0C12][D2C1C12][D2C1S12][D2C2C12][D2C2S12] . . . 


[TH][0] 

-[th][uzhoi] 

[th][uzhoi] 

-[th][uZH02] 

[th][uZH02] 


Substituting equation (K4) into equations (K3) results in 

[m - [MtrtllW + [c]W 1 [B2][E2]]] {X2} - {P2} - MUrVl) 

{xl} = [a] _1 {p1} - [[a]~ 1 [Bl][El] + [a] -1 [b2][e2]]{x2} 


In this form of equations (K3) , the effects of structural changes in the airframe 
contained solely in [d], [ei], and [e2]. 


(K5b) 

(K6a) 

(K6b) 

are 
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ft 


COMMENTS ON ORDER AND BANDEDNESS OF REDUCED HARMONIC BALANCE EQUATIONS 

Original Equations 

The indefinite form of the original reduced harmonic balance equations is given 
in equation (65) . As has been previously discussed, to effect a harmonic balance 
solution in practice, the series in equation (48) must be truncated. The assumption 
of a definite (i.e., finite) form for the series given by equation (48) results in a 
definite form of equation (65); that is, the coefficient matrix becomes a matrix of 
finite order. Truncations in the series represented by equations (2a) to (2c) do not 
affect the order of the coefficient matrix of equation (65) , although truncations in 
these series can affect the accuracy of the computations of the elements of the coef- 
ficient matrix. On the other hand, truncations in the series given in equations (2a) 
to (2c) do affect the bandwidth of the coefficient matrix, a characteristic which can 
be an important consideration relating to the computational effort required to solve 
equation (65) . Explicit formulas can be given for maximum values of order and band- 
width of the coefficient matrix. Let NT denote the number of terms retained in the 
series of equation (48). (For the example given by eq. (57), NT = 9.) Then the 
maximum order of the coefficient matrix in equation (65) is NT (NVB + NZH) by 
NT (NVB + NZH). Let MH denote the maximum harmonic number appearing in the series 
given by equations (2a) to (2c) . Then the resulting coefficient matrix has the 
banded form indicated in figure LI with maximum bandwidth given by 
2[l + 2 (MH) ] (NVB + NZH). 


Equations Rearranged for Solution by Matrix Partitioning 

Consistent with previous discussion related to maximum order and bandwidth of 
definite forms of equation (65) , the maximum order of the coefficient matrix in equa- 
tion (K2) is NT(NVB + NZH) by NT(NVB + NZH). However, because of the row and column 
interchanges necessary to transform equation (65) to the form given in equation (Kl) , 
the coefficient matrix of equation (K2) has the banded form depicted in figure L2 . 

The maximum bandwidths of each of the four submatrices of the coefficient matrix in 
equation (K2) are indicated. 
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ADDITIONAL COMPUTING SEQUENCES 

Two additional computing sequences for calculating airframe vibrations are 
described here to supplement the basic computing sequence given in the main text. 

The first sequence includes steps for reanalysis by the matrix partitioning method of 
appendix E. The second is a sequence which is appropriate when the rotor is repre- 
sented by impedances as described in appendix H. For each sequence, a block diagram 
is presented in which the tasks to be performed are arranged in three columns to dis- 
tinguish work done by rotor analysts, work done by airframe analysts, and joint work. 


Computing Sequence With Reanalysis 

Figure Ml shows a block diagram indicating the sequence of tasks for calculating 
airframe vibrations including reanalysis when members are varied. Supplementary notes 
follow, keyed by number to the individual blocks: 

1. Specify the linearized rotor model : The specification is made by providing 

the coefficient matrices in the Fourier series expansions of [mr] , [cr] , [kr] , and 
{QRO} represented by equations (2) . 

2. Identify rotor-airframe connection points : The rotor and airframe analysts 

agree on an arrangement of discrete points at which to designate connections. The 
displacements and rotations at these points needed to express connectivity are identi- 
fied. This defines the vector {c} discussed in conjunction with equations (C2) . 

3. Partition the rotor equations to isolate interface variables : Among the rotor 

variables {dvr} appearing in the linearized rotor equations, a subset {dvh} is identi- 
fied as the variables which explicitly appear in expressions characterizing connec- 
tions of the rotor model to the airframe model (eq. (C2a) ) . The matrices specifying 
the linear rotor model are partitioned accordingly. See discussion of equation (3) 
and following. 

4. Form rotor contributions to the coupling equations : The matrices [thr] and 

[ger] are formed. See equations (C2a) and (Clla) and related discussion. 

5. Identify harmonics in assumed solution : The steady-state solution given by 

equation (48) is specialized to a definite form agreed on by the rotor and airframe 
analysts. See equation (57) for an example of a definite form. 

6. Specify loads impinging directly on airframe : The distribution of external 

oscillating loads impinging directly on the airframe must be described to enable the 
airframe analyst to define the vector {l} appearing in equation (5) . Static compo- 
nents may be included, but it is not necessary to include them. Loads transmitted 
from the rotor to the airframe through the rotor-airframe connections are not included 
in this specification. 

7. F orm original finite-element model : A finite-element code is employed to 

generate the matrices [ma] , [ca] , and [ka], appearing in equation (5). Also generated 
are the load coefficients {LOC}, {Lis}, {Lie}, ..., appearing in equation (6). 
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8 . Partition the airframe equations to isolate the interface variables : Among 

the airframe variables {z}, a subset {zh} is identified to be the variables which 
explicitly appear in the expressions characterizing connections of the airframe model 
to the rotor model (eq. (C2b) ) . The matrices characterizing the airframe model are 
partitioned accordingly. See discussion of equation (7) and following. 

9. Form airframe contribution to the coupling equations : The matrix [tHa] is 

formed. See equation (C2b) and related discussion. 

10. Form the coupling equations : Compute [Tl] from equation (Clla) . Compute 

[T2] and [T3] from equation (Cllb) . Compute [th] from equation (C14a) . Compute [tc] 
from equation (C14b) . This establishes the coupling equations (9a) and (9b) . 

11. Modify the coupling equations : Identify the submatrices [tcd] and [tci] 

associated with the matrix [tc]. See equation (10) and related discussion. Identify 
the submatrices [thd] and [thi] associated with the matrix [th] . See equation (12) 
and related discussion. Compute the matrix [th] from equation (14) . The coupling 
equations are thus reduced to a single equation, equation (13) . 

12. Modify the airframe equations : The matrix [ti] in equation (15) is formed 

from equation (11) , as discussed in the text. Operating with [ti] on the partitioned 
form of the original finite-element equations, as shown by equation (16) , establishes 
the modified form of the airframe equations (eq. (21) ) . Note that this task has been 
to some extent routinized in the NASTRAN code through the feature called multipoint 
constraints . 

13. Partition airframe equations to isolate members to be varied : Identify the 

airframe members which are candidates for variation. Partition equation (60) to 
isolate these members in accordance with the discussion associated with equation (El) . 

14. Carry out matrix operations which are unaffected by member change s: For the 
Nth harmonic (N ^ 0) use equation (E7b) to compute the right side of equation (E2a) 
and compute the term 


[A13] 

[A23] 

[B13] 

[B23] 


-[B13] 
-[B23] 
[A13] 
[A23 j 


[A33] 

-[B33] 

-1 

[A3l] 

[A32] 

-[B3l] 

-[B32] 

[B33] 

[A33] 


_[B3l] 

[B32] 

[A3l] 

[A32] 


which appears in the coefficient matrix of equation (E2a) . Note that intermediate 
steps in this process yield the terms 


[A33] 

-[B33] 

-1 

[A3l] 

[A32] 

-[B31] 

-[B32] 

[B33] 

[A33] 


[B3l] 

[B32] 

[A3l] 

[A32] 
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and 


[A33] 

-[B33] 

-1 

[F 3] 

[B33] 

[A33] 


[G3] 


which appear in equation (E2b) . For the zeroth harmonic, use equation (ElOb) to com- 
pute the right side of equation (E9a) and compute the term 


[A13] 


[A33] 1 [[A3l] 


[A23] 



which appears in the coefficient matrix of equation (E9a) . Note that intermediate 
steps in this process yield the terms 


[A33] 1 [j>3l] [A32]j 

and 

[a33] _1 [g3] 

which appear in equation (E9b) . 

15. Solve for [xi], [X2], [x3], [Yl], [y2], [y 3]: For the nonzero harmonics, 
equation (E2a) is solved for [xi], [X2], [Yl], and [y 2] and the result substituted 
into equation (E2b) to obtain [x3] and [y 3]. For the zeroth harmonic, equation (E9a) 
is solved for [Yl] and [Y2] and the result substituted into equation (E9b) to 
obtain [Y3] . 

16. Assemble harmonic forced responses : For the nonzero harmonics, extract the 

matrix 


[uzhin] 

[uzhon] 

[UZSIN2] 

[UZSON2] 

[UZSIN3] 

[UZSON3] 

-[uzhon] 

[uzhin] 

-[UZSON2] 

[UZSIN2] 

-[UZSON3] 

[UZSIN3] 
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£ 5 " 


from the matrix 


[XI] 

[X2] 

[X3] 

[VI] 

[Y2] 

[Y3] 


as indicated by equation (E7a) (the latter matrix was computed in task 15) . Appro- 
priate rearrangement of the rows of the extracted matrix (reversing the reordering of 
rows and columns employed in arriving at equation (El)) directly yields the matrix of 
harmonic forced responses shown in sketch C. For the zeroth harmonic, extract the 
matrix 

[UZHIO] 

[UZSI02] 

[UZSI03] 

from the matrix 


[VI] 

[Y2] 

[Y3] 


as indicated in equation (ElOa) (the latter matrix was computed in task 15) . Appro- 
priate rearrangement of the rows of the extracted matrix yields the matrix of 
responses shown in sketch D. 

17. Complete airframe trim solution : The nonzero portions of the airframe trim 

displacement vector and the trim force vector are given by equations (D6) and (D8) . 
The displacements 


{ZHFNS} 
{ ZHFNC } 


and 


J {zSFNS} 1 
UZSFNC} J 


which appear in these equations are computed by extracting the third column of equa- 
tion (E7a) (obtained in task 15) and rearranging rows. This result and the results of 
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task 16 are used to compute the trim displacements and forces given by equations (D 6 ) 
and (D 8 ) . The corresponding results for the zeroth harmonic are computed from equa- 
tions (D9c) and (D9d) . 

18. Form red u ced harmonic balance equations ; The matrices [DX 4 X 2 X 3 X 4 ] are com- 
puted using equation (52) and the algorithm of appendix G. Equations (67) and the 
algorithm of appendix J are used to compute the coefficient matrix of equation (65) . 
In this process, the only rows and columns generated are those corresponding to 
retained harmonics (see task 5) . 

19. Solve reduced harmonic balance equations : Equation (65) is solved. 

20. Recover airframe responses : Equations (63) and (71) are used. If interface 

forces are desired, use equations (72) and (73) . 

21. Vary members ; Changes are made to any of the members which were previously 
identified as candidates for variation (see task 13) . Go back to task 15. 


Computing Sequence When Rotor Is Represented by Impedances 

■C 

Figure M2 shows a block diagram indicating the sequence of tasks for calculating 
airframe vibrations when the rotor system is represented by impedances, as discussed 
in appendix H. Supplementary notes follow, keyed by number to the individual blocks: 

1. Identify harmonics in assumed solution : The steady-state solution given by 

equation (48) is specialized to a definite form agreed on by the rotor and airframe 
analysts. See equation (57) for an example of a definite form. 

2. Identify rotor-airframe connection points : The rotor and airframe analysts 

agree on an arrangement of discrete points at which to designate connections . The 
displacements and rotations at these points needed to express connectivity are identi- 
fied. This defines the vector {c} discussed in conjunction with equations (C2) . 

3. Identify interface variables in nonlinear rotor equations : Among the rotor 

variables {VR} appearing in the nonlinear rotor equations, a subset {vh} is identified 
as the variables which explicitly appear in the expressions characterizing connections 
of the rotor model to the airframe model (eq. (Cla) ) . 

4. Specify the linearized rotor model : The specification is made by providing 

the impedance matrix (the coefficient matrix) in equation (H5) and the vector of 
forces appearing on the right side of that equation. 

5. Form rotor contributions to the coupling equations : The matrices [thr] and 

[ger] are formed. See equations (C2a) and (Clla) and related discussion. 

6 . Specify loads impinging directly on airframe : The distribution of external 

oscillating loads impinging directly on the airframe must be described to enable the 
airframe analyst to define the vector {l} appearing in equation (5) . Static compo- 
nents may be included, but it is not necessary to include them. Loads transmitted 
from the rotor to the airframe through the rotor-airframe connections are not 
included in this specification. 
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7. Form original finite-element model ; A finite-element code is employed to 
generate the matrices [ma], [ca], and [ka], appearing in equation (5) . Also gener- 
ated are the load coefficients {LOC}, {Lis}, {Lie}, . .., appearing in equation (6). 

8. Partition the airframe equations to isolate the interface variables ; Among 
the airframe variables {z}, a subset {zh} is identified to be the variables which 
explicitly appear in the expressions characterizing connections of the airframe model 
to the rotor model (eq. (C2b) ) . The matrices characterizing the airframe model are 
partitioned accordingly. See discussion of equation (7) and following. 

9. Form airframe contribution to the coupling equations : The matrix [tha] is 

formed. See equation (C2b) and related discussion. 

10. Form the coupling equations : Compute [Tl] from equation (Clla) . Compute [T2] 

and [T3] from equation (Cllb) . Compute [th] from equation (C14a) . Compute [tc] from 
equation (C14b) . This establishes the coupling equations (9a) and (9b) . 

11. Modify the coupling equations ; Identify the submatrices [TCD] and [TCl] 
associated with the matrix [tc]. See equation (10) and related discussion. Identify 

the submatrices [thd] and [THl] associated with the matrix [th] . See equation (12) 
and related discussion. Compute the matrix [th] from equation (14) . The coupling 
equations are thus reduced to a single equation, equation (13) . 

12. Modify the airframe equations : The matrix [ti] in equation (15) is formed 

from equation (11) , as discussed in the text. Operating with [ti] on the partitioned 
form of the original finite-element equations, as shown by equation (16) , establishes 
the modified form of the airframe equations (eq. (21)). Note that this task has been 
to some extent routinized in the NASTRAN code through the feature called multipoint 
constraints. 

13. Compute harmonic forced responses : Equation (60) is solved for each nonzero 

harmonic of interest (see task 1) and the results are assembled as shown in sketch C. 
For any zeroth harmonic responses, equation (61) is solved and the results assembled 
as shown in sketch D. 

14. Compute airframe trim solution ; Equation (D3) is solved for { ZHFNs} , {ZSFNS}, 
{ZHFNC}, and {ZSFNC}. The vectors representing trim displacements and forces are com- 
puted from equations (D6) and (D8) . For the zeroth harmonic, equations (D9b) to (D9d) 
are used. 

15. Form reduced harmonic balance equations : Equation (H10) and the algorithm of 

appendix J are used to compute the coefficient matrix in equation (H9) . In this 
process, the only rows and columns generated are those corresponding to retained 
harmonics (see task 1) . 

16. Solve reduced harmonic balance equations : Equation (H9) is solved. 

17. Recover airframe responses : Equations (H8) and (71) are used. If interface 

forces are desired, use equations (72) and (73) . 
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ROTOR ANALYSIS 


JOINT ANALYSIS 


AIRFRAME ANALYSIS 



Figure Ml.- Block diagram indicating the sequence of tasks for calculating 


airframe vibrations including reanalysis when members are varied. 


Ill 
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ROTOR ANALYSIS JOINT ANALYSIS AIRFRAME ANALYSIS 



Figure M2.- Block diagram indicating the 
airframe vibrations when rotor is 


sequence of tasks for calculating 
represented by impedances. 
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SYMBOLS 


Symbols used throughout the main text and appendices in this paper are included 
in this list. However, for the sake of brevity and simplicity, symbols used and 
defined in one section of the main text or in one appendix are not included. The 
orders (or maximum orders) of all matrices appearing in this paper are given in 
table I after the symbol list. 


The symbols OC, NC, and NS attached to a matrix symbol denote the following: 

OC constant (zeroth harmonic) term in Fourier series expansion 

NC Nth harmonic (N ^ 0) cosine term in Fourier series expansion 

NS Nth harmonic (N ^ 0) sine term in Fourier series expansion 

For example, {DVBOC}, {DVBNC}, and {DVBNS} are the constant term. Nth harmonic cosine 
term, and Nth harmonic sine term in the Fourier series expansion of {dVb} (see 
eq. (48) ) . 

"\ 

[A11],[A12],[A13], 

[A2l],[A22],[A23], \ 

[A3l],[A32],[A33] 

J 

[B11],[B12],[B13], 

[B2l],[B22],[B23], \ 

[B3l],[B32],[B33] 

J 

[5] matrix characterizing linear equations of constraint (see eq. (B8) ) 

[ca] viscous damping matrix for modified airframe finite-element model 

(see eq. (19) ) 


submatrices appearing in equation (El) 


submatrices appearing in equation (El) 


[CAH], [CA12],') 
[CA21], [CA22] 




submatrices of [ca] (see eq. (21)) 


[ca] viscous damping matrix for original airframe finite-element model 

(see eq. (B3) ) 


[CA11],[CA12], 

[CA21],[CA22] 


submatrices of [ca] (see eq. (B6) ) 


[cr] matrix appearing in linearized generalized force expression for rotor 

(see eq. (A19) ) 
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[CR11],[CR12] 
[CR21], [CR22] 


submatrices of [cr] (see eq. (A25)) 


{c} vector containing displacements and rotations associated with points con- 

necting rotor model to original airframe finite-element model 

{dvb} vector of small deviations from trim state {VBO} 

{dvh} vector of small deviations from trim state {VHO} 

{DVR} vector of small deviations from trim state {vro} 


[dx^X^xJ 


[0x^2X3X411] 

[dx 1 X 2 X 3 X 4 12] 

[DXjX^X^l] 

[0X3X3X3X422] 


matrix defined by equation (52) 


submatrices of [DX^X^X^X^] ( see e< 3- (51) ) 


{dz} vector of small deviations from trim state {zo} 

{DZH} vector of small deviations from trim state {zho} 

{dZS} vector of small deviations from trim state {zso} 

[eX 3 X 2 X 3 X 4 ] designation for typical submatrix appearing in coefficient matrix 

of equation (65) 

{f} vector of generalized force expressions for modified airframe finite-element 

model (see eq. (19) ) 

{f} vector of generalized force expressions for original airframe finite-element 

model (see eq. (B3) ) 

[Fl],[F2],[F3] submatrices appearing in equation (El) 

{fh} vector of generalized force expressions for modified airframe finite-element 

model associated with {zh} 

{fho} trim solution value for {fh} 

{fhtr} total interface forces required to suppress mean interface displacements 
(see eq. (56) ) 

{fo} trim solution for {f} 

{fs} vector of generalized force expressions for modified airframe finite-element 

model associated with {zs} 
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{FSO} trim solution value for {fs} 

[FUIN], [fuOn] matrices of in-phase and out-of-phase forces associated with 

Nth harmonic (N ^ 0) forced response (see eq. (60) ) 

[fuio] matrix of in-phase forces associated with zeroth harmonic forced response 
(see eq. (61) ) 

[FX 1 X 2 X 3 X^] designation for typical submatrix appearing in coefficient matrix in 

equation (72) and defined in equation (73) 

constant (zeroth harmonic) term in Fourier series expansion of the resultant 
force imposed on the airframe by the rotor 

Nth harmonic (N ^ 0) cosine term in Fourier series expansion of the 
resultant force imposed on the airframe by the rotor 

Nth harmonic (N ^ 0) sine term in Fourier series expansion of the resultant 
force imposed on the airframe by the rotor 

[G1],[G2],[G3] submatrices appearing in equation (El) 

[GER] matrix reflecting resultant of all elementary row operations 

H,L,P,Q indices appearing in equation (52) 

[ka] stiffness matrix for modified airframe finite-element model (see eq. (19) ) 

[KAll], [KA12] 1 

> submatrices of [ka] (see eq. (21)) 

[KA2l],[KA22 ] J 

[ka] stiffness matrix for original airframe finite-element model (see eq. (B3) ) 


submatrices of [ka] (see eq. (B6)) 


appearing in linearized generalized force expression for rotor 
eq. (A19) ) 


submatrices of [kr] (see eq. (A25)) 


of loads directly applied to modified airframe finite-element model 
sources which are external to both the airframe and the rotor 
mechanical system (see eq. (19)) 

{LH} vector of modified external loads associated with {zh} 


[KAll], [KA12] 
[KA21], [KA22] 

[kr] matrix 

(see 

[KRll], [KR12] 
[KR2l],[KR22] 

{l} vector 

from 




{fzhoc} 

{fzhnc} 

{fzhns} 
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{LS} vector of modified external loads associated with {zs} 

{£} vector of loads directly applied to original airframe finite-element model 

from sources which are external to both the airframe and the rotor 
mechanical system (see eq. (B3) ) 

{lh} vector of original external loads associated with {zh} 

{ls} vector of original external loads associated with {zs} 

[ma] mass matrix for modified airframe finite-element model (see eq. (19)) 


[MAll], [MA12]! 


[MA2l], [MA22] 


7 


submatrices of [ma] (see eq. (21) ) 


[ma] 


mass matrix for original airframe finite-element model (see eq. (B3) ) 


[MA11],[MA12] 1 
[ma21], [MA22] j 


submatrices of [ma] (see eq. (B6) ) 


[mr] matrix appearing in linearized generalized force expression for rotor 

(see eq. (A19) ) 


[MRU] [MR12] 
[MR2l] [MR22] 


submatrices of [mr] (see eq. (A25)) 


[MROCll] constant (zeroth harmonic) term in Fourier series expansion of [MRll] 
[MR1S11] first harmonic sine term in Fourier series expansion of [MRll] 
[MRlCll] first harmonic cosine term in Fourier series expansion of [MRll] 

N designates general Nth harmonic (N ^ 0) term 

NI number of elements in (c) 

NRC number of rows in [c] 

NT number of terms retained in equation (48) 

NV number of elements in {v} 

NVB number of elements in {vb} 

NVC number of elements in {vc} 

NVH number of elements in {vh} 
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NVR 

number 

of 

elements 

in 

{vr} 

NZ 

number 

of 

elements 

in 

{z} 

NZ 

number 

of 

elements 

in 

{z} 

NZH 

number 

of 

elements 

in 

{zh} 

NZH 

number 

of 

elements 

in 

{zh} 

NZS 

number 

of 

elements 

in 

{zs} 

NZS 

number 

of 

elements 

in 

{zs} 

NZSF 

number 

of 

variables 

. in {zs} corresponding to masses and structural members 


which are not to be varied 

NZSV number of variables in {zs} corresponding to masses and structural members 

which are to be varied 

{q} vector of all generalized force expressions for rotor system 

Cqb} vector of rotor generalized force expressions associated with {vb} 

{qbo} trim solution value for {qb} 

{qc} vector of rotor generalized force expressions associated with {vc} 

{QCO} trim solution value for {qc} 

{qh} vector of rotor generalized force expressions associated with {vh} 

{QHTR} specified mean value of {qh} 

{qr} vector of rotor generalized force expressions associated with {vr} 

{QRO} trim solution value for {QR} 

{QZHOC} constant (zeroth harmonic) term in Fourier series expansion of new airframe 
variables (see eq. (62a) ) 

{QZHNC} Nth harmonic (N ^ 0) cosine term in Fourier series expansion of new airframe 
variables (see eq. (62b) ) 

{qzHNs} Nth harmonic (N ^ 0) sine term in Fourier series expansion of new airframe 
variables (see eq. (62b) ) 

[tc] matrix characterizing linear equations of constraint on original airframe 

finite-element model inherent in formulation of rotor model (see 
eq. (C13b) ) 

[tCD],[tci] submatrices of [tc] consistent with the partitioning of {zh} into 
{ZHD} and {ZHl} 
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[th] matrix characterizing the linear relation connecting {DVH} to {ZHl} in the 

modified coupling equations (see eq. (13)) 

[th] matrix characterizing the linear relationship connecting {DVH} to {zh} in 

the original coupling equations (see eq- (C13a) ) 

[tha] matrix characterizing the linearized relation connecting {c} to {zh} 

[THR] matrix characterizing the linearized relation connecting {c} to {DVH} 

[thd] , [THl] submatrices of [th] consistent with partitioning {zh} into {zhd} 

and (ZHl) 

[Tl] matrix defined by equation (15) 

[Tl] , [T2] , [T3] submatrices appearing in equation (C9) 
t time 

submatrices of matrix of harmonic forced responses for Nth harmonic 
(N ^ 0) (see eq. (60) and sketch C) 


submatrices of matrix of harmonic forced responses for zeroth 
harmonic (see eq- (61) and sketch D) 



[uzhin] 

[uzhon] 

[uzsin] 

[uzson] 



parameters appearing in equation (52) and computed by the algorithm 
in appendix G 


{VB} vector of rotor generalized coordinates not contained in {vh}, {vc}, and VT 

{VBO} trim solution value for {vb} 

{vc} vector of rotor generalized coordinates which describe actions of the flight 

control system 

{vco} trim solution value for {vc} 

{vh} vector of rotor generalized coordinates which explicitly appear in expres- 

sions characterizing mechanical connections of the rotor system to the 
airframe system 

{VHO} trim solution value for {vh} 
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{vr} vector representing {Vb} and {vh} collectively 

{VRO} trim solution value for {vr} 


VT rotor generalized coordinate representing shaft rotation angle 

{v} vector containing all generalized coordinates of rotor mechanical system 

6w virtual work for airframe system 

ci 

6w r virtual work for rotor system 

6w virtual work of airframe system and rotor system combined 


[X1],[X2],[X3] 


submatrices appearing in equation (El) 


Xj,Yj,Zj inertial axis system 

X 1 X 2 X 3 X 4 four-character indexing system: the first character is an integer 

indicating row harmonic, the second character is the letter S or 
C indicating sine or cosine, the third character is an integer 
indicating column harmonic, and the fourth character is the 
letter S or C indicating sine or cosine 


[Y1],[Y2],[Y3] 


submatrices appearing in equation (El) 


{z} vector of all airframe generalized coordinates for modified finite-element 

model (see eq. (19)) 

{z} vector of all airframe generalized coordinates for original finite-element 

model (see eq. (B3) ) 

{zh} vector of airframe generalized coordinates which explicitly appear in 

expressions characterizing mechanical connections of the modified airframe 
finite-element model to the rotor model 

{ZH} vector of airframe generalized coordinates which explicitly appear in 

expressions characterizing mechanical connections of the original airframe 
finite-element model to the rotor model 

{zhd} variables in {zh} which are rendered dependent 

{zHl} variables in {zh} which remain independent 

{ZHO} trim solution value for {zh} 

{zo} trim solution for {z} 

{zs} vector of all airframe generalized coordinates of modified finite-element 

model not contained in {zh} 

{zs} vector of all airframe generalized coordinates of original finite-element 

model not contained in {zh} 
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{zso} 


trim solution value for {zs} 


{ZHFNC}^ 

{ZSFNC} 

{ZHFNS} 

{ZSFNS}^ 

{ZHFOC} 
{ZSFOC} 

a, 6/Y 



6 ( ) 

[ ] 

{ } 

[ ] T 
[ r 1 

[i] 

[ 0 ] 


Nth harmonic (N ^ 0) cosine and sine responses of airframe due to Nth 
harmonic external loads acting on the airframe (see eq. (D3) ) 


static (zeroth harmonic) response of airframe due to static external loads 
acting on the airframe (see eq. (D9b) ) 

parameters appearing in equations (67c) , (67d) , and (H10) computed by the 

algorithm given in appendix J 

virtual variation 

rectangular matrix 

column matrix or vector 

matrix transpose 

matrix inverse 

unit matrix 

null matrix 


(") 

C) 

C) 


second time derivative, 
first time derivative. 


d 


2 


dt 


2 


d_ 

dt 


airframe term associated with original finite-element model 
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TABLE I.- ORDER OF MATRICES 



Order 


Order 

Matrix 



Matrix 



Rows 

Columns 

Rows 

Columns 



1 1 

> 

1 I 

PJ 

(NVB) (NT) 

(NVB) (NT) 

[c] a 

(NZH) (NT) 

(NVB) (NT) 

[All] 

NZH 

NZH 

[ 5 ] 

NRC 

NZ 

[A12] 

NZH 

NZSV 

[ca] 

NZ 

NZ 

[A13] 

NZH 

NZSF 

[CAll] 

NZH 

NZH 

[A21] 

NZSV 

NZH 

[CA12] 

NZH 

NZS 

[A22] 

NZSV 

NZSV 

[CA2l] 

NZS 

NZH 

[A23] 

NZSV 

NZSF 

[CA22] 

NZS 

NZS 

[A3l] 

NZSF 

NZH 

[cl] 

NZ 

NZ 

[A32] 

NZSF 

NZSV 

[CAll] 

NZH 

NZH 

[A33] 

NZSF 

NZSF 

[CA12] 

NZH 

NZS 




[CA2l] 

nzs 

NZH 

[B] 3 

(NVB) (NT) 

(NZH) (NT) 

[CA22] 

NZS 

NZS 

[Bl] 3 

(NVB) (NT) 

(NVH) (NT) 

[cr] 

NVR 

NVR 

[B2] a 

(NVB) (NT) 

(NVH) (NT) 

[CRll] 

NVB 

NVB 

[Bll] 

NZH 

NZH 

[CR12] 

NVB 

NVH 

[B12] 

NZH 

NZSV 

[cR2l] 

NVH 

NVB 

[B13] 

NZH 

NZSF 

[CR22] 

NVH 

NVH 

[B2l] 

NZSV 

NZH 

rcRoci 

NVR 

NVR 

[B22] 

NZSV 

NZSV 

Tcrnc] 

NVR 

NVR 

[B23] 

NZSV 

NZSF 

Tcrns] 

NVR 

NVR 

[B3l] 

NZSF 

NZH 

{c} 

NI 

1 

[B32] 

NZSF 

NZSV 




[B33] 

NZSF 

NZSF 





a Order indicated for matrix is maximum possible value. 


TABLE X . - Continued 



Order 


Order 

Matrix 



Matrix 



Rows 

Columns 

Rows 

Columns 



CD ] 3 

(NZH) (NT) 

(NZH) (NT) 

{dzsnc} 

NZS 

1 

{dvb} 

NVB 

1 

{dzsns} 

NZS 

1 

{dvboc} 

NVB 

1 

[El] 3 

(NVH) (NT) 

(NZH) (NT) 

{dvbns} 

NVB 

1 

[E 2 ] a 

(NVH) (NT) 

(NZH) (NT) 

{dvbnc} 

NVB 

1 

[EXiX 2 x 3 x 4 ] 

NVB + NZH 

NVB + NZH 

{dvh} 

NVH 

1 

[ex-,x 9 x-,x a ii] 

NVB 

NVB 

{dvhoc} 

{dvhns} 

NVH 

1 

[EX4X2X3X412] 

NVB 

NZH 

NVH 

1 

NZH 

NVB 


[eX-,X 9 X,X^ 21 ] 

{dvhnc} 

{dvr} 

NVH 

1 

[EX-,X ? X ? X d 22 ] 

NZH 

NZH 

NVR 

1 




[DX4X2X3X4] 

NVR 

NVR 

{f} 

NZ 

1 

[0X2X2X3X411] 

NVB 

NVB 

{f} 



NZ 

1 

[0X4X2X3X412] 

NVB 

NVH 

[FI ] 3 

NZH 

2 (NZH) + 1 

[0X4X3X3X421] 

NVH 

NVB 

[F2] 3 

NZSV 

2 (NZH) + 1 

[0X4X3X3X422] 

NVH 

NVH 

[F 3 ] a 

NZSF 

2 (NZH) + 1 

{dz} 

NZ 

1 

{fdz} 

NZ 

1 

{dzh} 

NZH 

1 

{fdzh} 

NZH 

1 

{dzhoc} 

NZH 

1 

{fdzhoc} 

NZH 

1 

{dzhns} 

NZH 

1 

{fdzhnc} 

NZH 

1 

{dzhnc} 

NZH 

1 

{fdzhns} 

NZH 

1 

{dzs} 

NZS 

1 

{fdzs} 

NZS 

1 

{dzsoc} 

NZS 

1 

{fdzsoc} 

NZS 

1 


a Order indicated for matrix is maximum possible value. 
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TABLE I.- Continued 


Order j 

Rows 

Columns 

NZS 

1 

NZS 

1 

NZH 

1 

NZH 

1 

NZH 

1 

NZH 

1 

NZH 

1 

NZH 

1 

NZ 

1 

NZS 

1 

NZS 

1 

NZH 

NZH 

NZH 

NZH 

NZH 

NZH 

NZH 

NZH 

NZH 

1 

NZH 

1 

NZH 

1 

NZH 

NZH 

NI 

1 

NI 

1 


Matrix 

[Gl] a 

[G2] a 

[G3] a 

[ger] 

[gx 1 x 2 x 3 x 4 ] 

[g] 

[hio] 

[hin] 

[hon] 

[ka] 

[kaii] 

[KA12] 

[KA2l] 

[KA22] 

[ka] 

[kaii] 

[KA12] 

[KA2l] 

[KA22] 

[kr] 


Order 

Rows 

Columns 

NZH 

2 (NZH) + 1 

NZSV 

2 (NZH) + 1 

NZSF 

2 (NZH) + 1 

NI 

NI 

NZH 

NZH 

NZH 

NZH 

NZ 

NZ 

NZ 

NZ 

NZ 

NZ 

NZ 

NZ 

NZH 

NZH 

NZH 

NZS 

NZS 

NZH 

NZS 

NZS 

NZ 

NZ 

NZH 

NZH 

NZH 

NZS 

NZS 

NZH 

NZS 

NZS 

NVR 

NVR 


Order indicated for matrix is maximum possible value 


TABLE I.- Continued 


Order 


Matrix 

Rows 

[krii] 

NVB 

[KR12] 

NVB 

[KR2l] 

NVH 

[KR22] 

NVH 

[kroc] 

NVR 

[krnc] 

NVR 

[krns] 

NVR 

x 1 x 2 x 3 x 4 ] 

NVH 


NZ 

NZ 

NZ 

NZ 

NZH 

NZH 

NZH 

NZH 

NZS 

NZS 

NZSV 

NZSF 

NZS 

NZSV 

NZSF 


Matrix 


{LSNS} 
{LSNS2} 
{ LSNS 3} 
{£} 
{loc} 
{lns} 
(lnc) 
{lh} 
{lhoc} 


{lhnc} 

{lhns} 

{ls} 

{lsoc} 


Order 


Rows Columns 


NZS 

NZSV 

NZSF 


NZ 


NZ 


NZ 


NZ 


NZH 


NZH 


NZH 


NZH 


NZS 


NZS 


NZS 


NZS 



























TABLE I.- Continued 



Order 


Order 

Matrix 

Rows 

Columns 

Matrix 

Rows 

Columns 

[MA2l] 

NZS 

NZH 

{qbons} 

NVB 

1 

[MA22] 

NZS 

NZS 

{Qc} 

NVC 

1 

[mr] 

NVR 

NVR 

{QCO} 

NVC 

1 

[mroc] 

NVR 

NVR 

{Qh} 

NVH 

1 

[mrnc] 

NVR 

NVR 

{qhoc } 

NVH 

1 

[mrns] 

NVR 

NVR 

{qhnc } 

NVH 

1 

[MRll] 

NVB 

NVB 

{QHNS} 

NVH 

1 

[MR12] 

NVB 

NVH 

{QHO} 

NVH 

1 

[MR2l] 

NVH 

NVB 

{QHOOC} 

NVH 

1 

[MR22] 

NVH 

NVH 

{qhonc} 

NVH 

1 

[MR0C11] 

NVB 

NVB 

{qhons} 

NVH 

1 

[mrncii] 

NVB 

NVB 

{ qhtr} 

NVH 

1 

[mrnsii] 

NVB 

NVB 

{QO} 

NV 

1 




{QR} 

NVR 

1 

{Pl} a 

(NVB) (NT) 

1 

{QRO} 

NVR 

1 

{P2} a 

(NZH) (NT) 

1 

{QROOC} 

NVR 

1 

{Q} 

NV 

1 

{QRONC} 

NVR 

1 

{QB} 

NVB 

1 

{ QRONS} 

NVR 

1 

{qbo} 

NVB 

1 

{QZHOC} 

NZH 

1 

{qbooc} 

NVB 

1 

{ qzhnc} 

NZH 

1 

{qbonc} 

NVB 

1 

{qzhns} 

NZH 

1 


a Order indicated for matrix is maximum possible value. 




TABLE I.- Continued 


Matrix 

Order 

Matrix 

Order 

Rows 

Columns 


Rows 

Columns 

[tc] 

NI 

- NVH 

NZH 


[UZSON3] 

NZSF 

NZH 


[tcd] 

NI 

- NVH 

NI - NVH 


[UZSI02] 

NZSV 

NZH 


[TCI] 

NI 

- NVH 

NZH + NVH - 

NI 

[UZSI03] 

NZSF 

NZH 


[th] 


NVH 

NVH + NZH - 

NI 

{vb} 

NVB 

1 


[tha] 


NI 

NZH 


{VBO} 

NVB 

1 


[thr] 


NI 

NVH 


(vc> 

NVC 

1 


[th] 


NVH 

NZH 


(vco} 

NVC 

1 


[thd] 


NVH 

NI - NVH 


(vh) 

NVH 

1 


[thi] 


NVH 

NZH + NVH - 

NI 

(VHO) 

NVH 

1 


[ti] 


NZH 

NZH + NVH - 

NI 

(vo> 

NV 

1 


[Tl] 


NVH 

NVH 


{vr} 

NVR 

1 


[T2] 


NVH 

NZH 


(VRO) 

NVR 

1 


1 — 1 

H3 

CO 

1 1 

NI 

- NVH 

NZH 


(v) 

NV 

1 


[uzhin] 


NZH 

NZH 


(wc> 

NZH 

1 


[uzhon] 


NZH 

NZH 


{ws} 

NZH 

1 


[uzsin] 


NZS 

NZH 






[uzson] 


NZS 

NZH 


[xi] 

NZH 

2 (NZH) + 

1 

[UZHIO] 


NZH 

NZH 


[X2] 

NZSV 

2 (NZH) + 

1 

[uzsio] 


NZS 

NZH 


[X3] 

NZSF 

2 (NZH) + 

1 

[UZSIN2] 


NZSV 

NZH 


{xi} a 

(NVB) (NT) 

1 


[UZSIN3 ] 


NZSF 

NZH 


(X2 } a 

(NZH) (NT) 

1 


[UZSON2] 


NZSV 

NZH 







a Order indicated for matrix is maximum possible value. 
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TABLE I.- Concluded 


Order 



Rows 

Columns 

a 

NZH 

2 (NZH) + 1 

a 

NZSV 

2 (NZH) + 1 

a 

NZSF 

2 (NZH) + 1 


NZ 

1 


NZ 

1 


NZH 

1 


NZH 

1 


NI - NVH 
NZH + NVH - NI 
NZH 
NZH 
NZH 
NZH 
NZ 



Order 


Rows 

Columns 

{zs} 

NZS 

1 

{zs} 

NZS 

1 

{ ZSFNS } 

NZS 

1 

{zSFNC} 

NZS 

1 

{ZSFOC} 

NZS 

1 

{ZSFNS 2} 

NZSV 

1 

{ZSFNS3} 

NZSF 

1 

{ ZSFNC 2 } 

NZSV 

1 

{ ZSFNC3 } 

NZSF 

1 

{ZSF0C2> 

NZSV 

1 

{ZSF0C3} 

NZSF 

1 

{zso} 

NZS 

1 

{zsooc} 

NZS 

1 

{ ZSONC} 

NZS 

1 

(ZSONS } 

NZS 

1 


Order indicated for matrix is maximum possible value. 
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